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1 Introduction 

La découverte en 1965 par Penzias et Wilson du fond de rayonnement cos- 
mologique diffus à 3 K a définitivement fait du modèle du Big Bang chaud 
le cadre général de la cosmologie moderne. Ces trois dernières décennies ont 
été marquées par de nombreux progrès tant théoriques qu'observationnels qui 
n'ont fait que confirmer sa validité, et de préciser petit à petit les détails de ce 
scénario. 

Parmi les avancées notables de ces dernières années, citons pour mémoire: 

• La réalisation de grands catalogues de galaxies. On a maintenant des 
catalogues tri-dimensionnels qui contiennent des milliers d'objets. Ces 
catalogues permettent de faire une véritable cosmographie de l'univers 
local. De nouveaux moyens sont apparus récemment pour compléter ce 
type de catalogues, la détermination des champs de vitesse cosmiques, 
permet ainsi d'accéder à des informations cinématiques très précieuses. 
Et depuis la fin des années 80, un nouveau moyen d'investigation est en 
train d'émerger, il s'agit des cartes de distorsion gravitationnelle. Elles 
permettent de visualiser les lignes de potentiel de la masse projetée. 

• le modèle de Matière Noire Froide (CDM pour Cold Dark Matter) qui 
est apparu au début des années 80 et qui a servi de point de référence 
(à défaut de devenir un modèle standard) pour tous les travaux sur le 
problème de la formation des grandes structures. 

• Le développement des théories inflationnaires. C'est encore un terrain très 
spéculatif, mais c'est le lieu très excitant où les concepts de la physique 
des hautes énergies rencontrent des exigences observationnelles de plus en 
plus fiables. 
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• Enfin, la détection en 1992 des fluctuations de température du fond de ray- 
onnement cosmologique par l'expérience satellitaire COBE/DMR (Smoot 
et al. 1992) a marqué un tournant pour la cosmologie: pour la première 
fois on avait une preuve directe de l'origine des grandes structures de 
l'univers. 

Dans ce cours je vais principalement m'intéresser au problème de la forma- 
tion des grandes structures. 

Quelques ouvrages de référence, 

• Relativité Générale: Weinberg, Gravitation and Cosmology, 1972; Landau 
et Lifschitz, Classical Theory of Fields, 1975 

• Univers primordial: Kolb et Turner The Early Universe, 1990, mais il est 
peu précis; revue de Brandanberger, Inflation and Cosmic Strings: two 
Mechanisms for producing Structure in the Universe, Int. J. Mod. Phys. 
A2 : 77, 1987. 

• Inflation: Linde, Particle Physics and Inflationary Cosmology, 1990; Lid- 
dle et Lyth, The Cold Dark Matter density perturbation, 1993, Physics 
Reports, 231, 1; Lidscy et al. Reconstructing the inflaton potential - an 
overview, 1997, Reviews of Modem Physics, 69, 2 

• Formation des grandes structures: Peebles The Large Scale Structure of 
the Universe, 1980; Principle of Physical Cosmology, 1993; V. Sahni & 
P. Coles, Approximation Methods for Non-linear Gravitational Clustering, 
1995, Physics Reports, 262, 1 

2 Eléments de cosmologie, l'univers homogène 

2.1 Les principes cosmologiques et cosmographie 

Le principe fondamental sur lequel repose la cosmologie, est l'hypothèse que 
l'univers est, au moins à grande échelle, homogène et isotrope. Cette idée est 
née au début de ce siècle. Elle est en fait, à l'origine, davantage motivée par 
des considérations philosophiques que des raisons physiques observationnelles. 
Einstein a été un des premiers, dès la naissance de le Relativité Générale, à tenter 
de construire une métrique susceptible de décrire l'ensemble de l'univers. L'idée 
générale de départ était que l'univers était non seulement homogène, isotrope 
mais aussi immuable. L'abandon de cette dernière hypothèse est d'abord due à 
la mise en évidence en 1929, par un travail acharné de E. Hubble, du mouvement 
de récession des galaxies. Cela dit, il faudra quand même attendre 1965, avec la 
découverte du rayonnement fossile à 3K par Penzias & Wilson pour que l'idée 
que l'univers ait traversé une phase chaude, donc le modèle du Hot Big Bang, 
s'impose face aux tenants d'un Univers immuable qui défendaient le modèle 
'steady state' (modèle stationnaire). 
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On a aujourd'hui des évidences assez solides en faveur de l'homogénéité et 
de l'isotropie. Pendant longtemps les relevés locaux d'étoiles ou de galaxies ont 
montré de fortes hétérogénéités. Les plus grands catalogues de galaxies arrivent 
enfin à mettre en évidence une homogénéité à grande échelle. Le Las-Campanas 
Redshift Survey qui contient plusieurs milliers d'objets l'illustre assez bien. 

L'évidence la plus forte reste sans doute l'isotropie du fond de rayonnement 
à 3K. Les écarts à l'isotropie ont maintenant été détectés. Le mouvement propre 
de notre galaxie induit un écart de l'ordre de 1CP 3 , mais les anisotropies propres 
sur cette surface de dernière diffusion ne sont que de l'ordre de 10 -5 . 

Tout cela contribue à valider le modèle d'un Univers homogène et isotrope 
à grande échelle. 

Dans le cadre de la Relativité Générale le principe cosmologique détermine 
une forme générale de la métrique. On peut se référer pour ce calcul au livre de 
Weinberg (1972). La forme générale de la métrique qui répond à ces conditions 
est donnée par la métrique dite de Robertson-Walker, 



ds 2 = -At 2 + a 2 {t) 



H-r 2 

9 +x 2 d6 2 +x 2 sin 2 #d0 2 



(1) 



où k est une constante. Si k est nul la partie spatiale de la métrique est plate. 
Sinon on a une courbure non nulle. 

On voit que cette métrique est indépendante de la position et qu'il n'y a 
aucune direction privilégiée. Mais cela ne préjuge pas a priori de la topologie 
globale de l'univers. Même quand k est nul ou négatif, on peut très bien avoir 
un Univers de taille finie en imaginant une topologie périodique (voir une revue 
récente de Lachièze-Rey et Luminet, gr-qc/9605010, Phys. Rept. 254 (1995) 
135). Notons que cela brise tout de même l'isotropie globale de l'Univers. De 
telles conditions périodiques n'ont pas de conséquences sur la dynamique, sur la 
résolution des équations d'Einstein, qui sont purement locales (c'est vrai aussi 
bien pour le facteur d'expansion, que pour la formation des structures, au moins 
tant que ces structures sont plus petites que la maille élémentaire). Par contre 
cette hypothèse a des conséquences observationnelles dès que l'on s'intéresse 
aux propriétés statistiques des champs cosmiques: la même portion se retrouve 
dans des directions différentes induisant de ce fait des propriétés statistiques 
spécifiques. Dans la suite de cet exposé je vais négliger une telle possibilité. 

Une conséquence du facteur d'expansion a, c'est la loi de Hubblc. En effet 
la vitesse relative de deux particules à la distance physique, r = a x, est donnée 
par 

d à . . . 

V = dt^ aX ' = a V + aX ' 
Cette relation montre que la vitesse moyenne des objets les uns par rapport 
aux autres est proportionnelle à la distance de ces objets. On appelle H cette 
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Figure 1: Le Las-Campanas Redshift Survey, Shectman et al. 1996 
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quantité, dite 'constante de HubbleQ, 



Il y a un terme supplémentaire qui est la vitesse particulière des particules. 
Cette vitesse est en moyenne nulle. Pour une galaxie elle est de l'ordre de 300 
km/s. 

Cette relation là est purement locale, si on veut faire des calculs plus précis 
il faut tenir compte de la courbure spatio-temporelle. Imaginons qu'on détecte 
un signal lumineux (d'une galaxie lointaine) à un instant t^. Il a été émis à 
l'instant t\ et à la distance x%. Pour relier t\ à x\ il faut intégrer le long de la 
géodésique suivie par les photons. Dans ce cas c'est très simple, 

^=0 = At--aHt)^- 2 . (4) 



(5) 



On a donc, 

''<!*_ f Xl dx 
ti W)~ Jo Vl-kx 2 ' 
Si la galaxie émettrice est au repos, x\ indépendant du temps, alors deux 
événements séparés par un intervalle de temps 8t\ au moment de l'émission 
sont séparés d'un intervalle de temps <5 to à la réception avec, 

o(*o)~o(ti)' W 

ce qui implique que les longueurs d'onde électromagnétiques Ai et Ao sont reliées 
par, 

Ao _ a(t ) , 
Ai a(hY [ ' 

Or dans un spectre d'objets lointains comme les galaxies, des raies spécifiques, 
d'absorption par exemple, se retrouvent assez facilement. L'expansion (donc 
éloignement relatif des objets) se traduit par un décalage vers le rouge de ces 
raies, (que je vais dans toute la suite appelé redshift). En mesurant le rap- 
port d'une longueur d'onde mesurée avec la longueur d'onde correspondante au 
repos, on mesure donc directement un rapport de facteurs d'expansion entre 
l'émission et la réception, 

Notons que si on interprétait ce décalage spectral comme une vitesse d'éloignement, 
cette vitesse serait de cz/(l + z). 

Le redshift est finalement la mesure de distance par excellence en cosmologie. 
Si on a tendance à utiliser les Mega-parsecs (un parsec vaut approximativement 

1 ce terme est bien mal approprié puisque cette quantité dépend en fait du temps 
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3.5 année- lumière) , on réintroduit finalement la constante de Hubble dedans et 
c'est bien de vitesses relatives que l'on parle. Pour cela on introduit la variable 
réduite, 

h= —, f . (9) 

100 kms^Mpc" 1 

Ainsi la distance typique des galaxies entre elles est de l'ordre de 5/i _1 Mpc, 
c'est à dire de 500 km/ s. La difficulté qu'on a à mesurer la constante de Hubble 
vient du fait qu'il faut trouver un moyen de calibrer une relation. Il faudrait 
une " règle" , et on en a pas de très bonne à notre disposition. 

La relation de proportionnalité entre la vitesse et la distance n'est vraie que 
dans l'Univers local. Dès que le redshift approche l'unité il faut préciser de 
quelle distance on parle. 

• La distance angulaire. Un objet de taille propre D donnée sera vu sous 
un angle S. Si cet angle est petit, alors on peut définir la distance d A par 

Ss f-. (10) 

Si on écrit ds 2 = pour la distance transversale d'un objet, donc à x\ 
fixe, on a dt = a{t\) x\{ti) S d'après la forme de la métrique. On a donc 
simplement 

d A {ti)=a{t 1 )x 1 (h). (11) 

• La distance lumineuse. Cela correspond au fait que physiquement la lumi- 
nosité apparente d'une source décroît avec la distance. Plus précisément si 
^abs. est la luminosité absolue d'un objet et L app sa luminosité apparente, 
alors on peut définir par, 

Lahs. — iapp. 47rd|, (12) 

et on trouve 



d L (h) = -^d A (t 1 ). (13) 



a 2 (t ) 
a 2 (h) 

La meilleure manière d'obtenir cette relation est d'écrire la conservation 
de l'énergie des photons émis. 



• La distance du mouvement propre. Si l'angle, AS, parcouru par un objet 
de vitesse perpendiculaire V± pendant le temps (de l'observateur) Ai , on 
définir la distance du par, 

AS=Y^. (14) 

«M 

On a 

d M (ti) = ^\d A (t 1 ). (15) 
a(ti) 
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• La distance parallaxe. Cela correspond à la mesure d'une distance par 
effet de parallaxe (changement de la position de l'observateur). C'est un 
peu le contraire du cas de la distance angulaire. On trouve 



d P (ti) = a(t ) 



Xl 



y/l-kxl 



(16) 



Ce serait intéressant, mais cette méthode ne permet absolument pas d'accéder 
à des distance cosmologiques! 

Un paramètre qu'il peut être utile de considérer est le paramètre de décélération 



5o, 



aa 



On peut alors faire le développement perturbatif à t = t , 



a(t) = a a 



1 + H (t - t ) - ^ q H 2 (t - h) 2 + 



qui donne une relation perturbative entre la distance x\ et le redshift z, 



Xl 



1 



a H a 



-\(l + q )z 2 + ... 



(17) 



(18) 



(19) 



La comparaison distanec-redshift permet donc en premier lieu de mesurer le 
paramètre de décélération q . Ce qu'on vient de faire n'est que de la cosmo- 
graphie à partir de la métrique de R.W.. Pour faire mieux il faut écrire les 
équations dynamiques. 



2.2 Évolution du facteur d'expansion 

Les équations d'Einstein s'écrivent, 



R 



-8tt G + \ g^ v . 



(20) 



Je garde ici la possibilité d'avoir une constante cosmologique non nulle. Une 
telle métrique impose que le tenseur énergie-impulsion du fluide cosmique est 
isotrope (pas de pression anisotrope). Autrement dit, peut s'écrire 



T 



( P 



V 



\ 



p ) 



(21) 



A partir de la conservation de l'énergie impulsion on peut écrire (composante 
(0,0)), 

3 dp _ d 



dt dt 



[a 3 (P+P)} 



(22) 
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qui peut se réécrire, 

±(pa 3 ) = -3pa 3 . (23) 
dt 

Pour simplifier les discussions à venir je vais supposer qu'il existe une équation 
d'état de la forme, 

p = wp, (24) 

qui relie la pression locale à la densité, w est supposé constant. Évidemment 
cette hypothèse, vraie pour chacun des fluides pris séparément, n'est pas vraie en 
général pour le fluide total mais c'est vrai une grande partie du temps (sauf pen- 
dant les périodes de transition). Quand la densité d'énergie du rayonnement (ou 
des particules ultra-relativistes) domine sur celle des particules non-relativistes, 
on a alors w — 1/3. Si c'est la matière qui domine alors on a w — (poussières 
sans pression). 

Notons quand même que la conservation d'énergie s'applique pour chaque 
fluide cosmique séparément (si ils ne sont effectivement couplés que par gravi- 
tation). Cela donne l'évolution de la densité de chacun des fluides en fonction 
du facteur d'expansion. 

p 1 ex a -4 pour le rayonnement (25) 
Pmat. oc aT 3 pour la matière. (26) 

Les équations d'Einstein donnent alors les équations d'évolution du facteur 
d'expansion a. 

3 a = -A-kG p (l + 3w) a + Aa; (27) 



• 2 1 
a = — 



8"7rG o , A(r 
— — pa - ka+ — 



(28) 



Si on interprète A comme un fluide cosmique, on a une densité associée de 
px = —p\ = X/(8ttG). Noter que, compte tenues des équations de conservation, 
l'équation ( p8| ) est l'intégrale de la première et k joue là le rôle d'une constante 
d'intégration. Cela permet de donner une interprétation plus dynamique à k. 
Plus précisément, la première équation correspond à l'application du théorème 
de Gauss (ou plutôt de BirkhofF dans le cadre de la relativité générale) à une 
sphère de densité p+A, alors, si A = 0, cette sphère peut croître jusqu'à atteindre 
un rayon infini uniquement si k < 0. Si k > 0, il y a un rayon d'extension 
maximale. 

Pour avoir la dépendance en temps du facteur d'expansion il faut finalement 
résoudre les équations complètes. 
Si k — et A = 0, on trouve 

a(t) ex t 2/3 , puisque à 2 - - (29) 
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Figure 2: Application du théorème de Birkhoff au facteur d'expansion 



pour un Univers dominé par la matière et 

ait) ex t 1/2 , puisque à 2 ~ \ (30) 

a 

pour un Univers dominé par le rayonnement. Pour un Univers dominé par la 
matière quand A = on a une solution explicite pour la croissance du facteur 
d'expansion. Si k = +1 (Q > 1) on aQ 

et si k = — 1 on passe à des fonctions hypergéométriques en faisant 9 — *■ hp. 

Je reviens tout de suite sur l'équation d'évolution de a, Eq. (p8|), et des 
différents termes qui la composent. Si on suppose a priori que a croit avec le 
temps, on voit que les termes en p, k puis A dominent successivement le terme 
de droite. 

Pour simplifier la discussion on introduit deux quantités réduites, 

3 H 2 

il = pi p c avec p c = — — . (33) 

et 

A -5^. (34) 

On a alors 

+ A = f + 4- (35) 



'Dans toute la suite les indices font référence à des quantités prises au temps présent 



a/a„ 




tH, 



o 



Figure 3: Évolution du facteur d'expansion dans les cas k = — 1, k — ou 
k = +1. 



qo 



A . 



(36) 



Si k = la somme des 2 doit faire 1 par construction. Remarquons aussi que 
ces deux quantités sont en générale dépendantes du temps, fi = 1 et A = 
étant le seul point statique (intéressant). On peut remarquer que le paramètre 
de décélération q s'identifie à, 

2 

Dans le plan (f2, A) on a, certes, de plus en plus d'indications sur le point 
qu'occuperait actuellement l'Univers mais il faut bien avouer qu'aucune mesure 
de f2 ou A n'est vraiment probante et indiscutable. On sait quand même que 
Çî doit être compris entre 0.2 et 1. Pour A, c'est encore moins précis, disons 
qu'il est de l'ordre de quelques unités. Si A est trop grand la distance lumineuse 
n'est plus monotone, et finalement on a plus de Big-Bang. Si A est trop petit 
(négatif) l'Univers est très jeune. 

Dans ce plan, il y a deux segments qui sont plus particulièrement envisagés 
et qui correspondent à des a priori théoriques différents, 

• Le cas A = 0; 

• Le cas fl + A = 1 donc k = 0; 

Sur la figure (|J), j'ai représenté les deux directions d'instabilité, et la dépendance 
de f2 et A avec le temps dans ces deux directions. 

Avec une solution explicite pour le comportement du facteur d'expansion, 
on peut expliciter la dépendance des distances avec le redshift, 



xi(z) 



■qo + (go -!)(-! + V'2 go z + l) 
H a ql (1 + 2) 



(37) 
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No Big Bang 




Figure 4: Le plan fi- A et évolution des paramètres cosmologiques réduit fi et A 
dans les cas particulier où A = ou k = 0. 



Cette expression est valable pour A = 0. Evidemment on retrouve le développement 
précédent qui est valable pratiquement jusqu'à z = 1. Pour Slo = 1, on a 
q = 1/2, on a, 

Xl = 4- ( 1 - (38) 



ff V Vï+^/ 

Un concept intéressant à considérer est l'horizon d'une particule. Si on a un 
observateur à un temps to, on peut se demander qu'elle elle la distance maximale 
peut se trouver un événement observable, donc 1"' horizon de la particule". Cet 
horizon, xh, est naturellement défini par, 



dx 



o VT^k 



x J t 



dt 



(39) 



Cette intégrale converge si le facteur d'expansion ne croît pas trop vite à l'origine. 
On a, 



Avec Ho to = 2/3 on a, 



2 / a 



X H = 



Ho \cio 



1/3 



3/2 



(40) 



(41) 
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La taille de notre propre horizon est donc de l'ordre de 6000 h~ 1 Mpc. Une 
quantité auquelle on peut s'intéresser c'est la taille apparente de l'horizon d'une 
particule sur la surface de dernière diffusion. Cette angle est donnée par, 

La taille de l'horizon sur le ciel est donc de l'ordre de 2deg: si le facteur 
d'expansion a évolué de manière standard depuis l'origine du Big-Bang, des 
régions du ciel séparées de plus de cette distance n'ont donc jamais été en con- 
tact causal. 



2.3 Contraintes actuelles sur les paramètres cosmologiques 

Dans ce paragraphe je vais essayer de faire le point sur la détermination des 
paramètres cosmologiques. 

• La constante de Hubble. La mesure de la constante de Hubble est une 
entreprise très difficile. La maîtrise des biais systématiques est en parti- 
culier très ardue. Le satellite Hipparcos a apporté quelque lumière à la 
première étape de l'échelle cosmique, en mesurant la distance par par- 
allaxe de céphéides de notre voisinage. Le Hubble Space télescope s'est 
évertué à rechercher des céphéides dans des galaxies lointaines (Virgo). 
Les valeurs les plus souvent admises sont de l'ordre de 60 à 70 km/s/Mpc. 
Globalement les résultats récents indiquent quand même une réduction du 
problème de l'âge de l'Univers. 

• Le paramètre de densité, fL Les estimations dynamiques donnent des 
valeurs de l'ordre de 0.2 (pour les amas) à l'unité (flots cosmiques à grande 
échelle, voir Dekel 1994). Je reviendrai sur ces méthodes qui sont basées 
sur des études de structures. 

• La constante cosmologique. En plus des contraintes mentionnées plus 
haut, on s'attend à ce que A ne dépasse pas l'unité, sinon le volume de 
l'Univers serait tel qu'on aurait énormément d'objets lointains (quasars) 
qui subiraient des effets de lentille. 

• le paramètre de décélération q . La relation distance-redshift, à grand 
redshift, est sensible au premier ordre au paramètre de décélération qg. 
L'idée est de se servir de Super-Novae, qui sont à des redshift connus, et 
de déterminer leur distance en s'appuyant sur leur courbe de lumière (voir 
figure). Pour l'instant les contraintes ne sont pas forcément très solides. 

Pour donner un aperçu des dernières tendances, globalement il semble 
qu'un modèle où 

fi = 0.3 (43) 

A = 0.7 (44) 
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avec une courbure nulle n'est pas une hypothèse ridicule et commence à 
être prise sérieusement en compte. 

2.4 Histoire thermique de l'Univers 

Je poursuis mes rappels sur le modèle du bing-bang chaud en faisant un bref de- 
scriptif de l'histoire thermique de l'univers. Dans les modèles qui nous intéressent 
le facteur d'expansion est évidemment une fonction croissante du temps, ce qui 
implique que la densité et éventuellement la pression croissent quand on remonte 
dans le passé. 

On peut identifier un certain nombre d'étapes clefs. On peut comparer la 
densité de rayonnement à la densité de matière (ou du moins ce que l'on en 
sait). 

Pour le rayonnement, c'est la densité du corps noir cosmologique, dont la 
densité d'énergie est donnée par la loi de Stephan, 

p 7 = aT 4 , (45) 

où a est la constante de Stephan. Notez que c'est là une quantité très bien 
connue! (voir figure montrant l'émission du corps noir). La densité de matière, 
c'est par définition 

Pmat. = QoPc- (46) 

Celle-ci est nettement moins bien connue, à cause de fi et de Hq. Il n'empêche, 
on est aujourd'hui largement dominé par la matière, mais quand on remonte 
dans le passé il y a un moment où p 7 = p ma t- Cette équivalence a lieu quand 

a = a cq . (47) 

avec 

aeq = 40 000O /i 2 ' (48) 
A ce moment là la température est à peu près de 100 000 K. 

C'est un moment crucial pour l'évolution du facteur d'expansion, mais ce qui 
est aussi intéressant c'est ce qui se passe pour les baryons quand la température 
et la pression montent. A une température de quelques milliers de degrés la 
'recombinaison' a eu lieu: au delà de cette température électrons et protons se 
découplent. Il reste de cette époque le fond de rayonnement à 3 K. 

Je vais illustrer ce qui se passe en calculant explicitement le redshift de la 
recombinaison. Cela va illustrer la méthode générale mise en œuvre pour faire 
ce genre de calculs. 

On suppose qu'on a un équilibre chimique entre les différentes espèces qui 
interviennent, 

e +p <-> H +-f. (49) 
On écrit l'équilibre des potentiels chimiques et on obtient, 
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z 



0.02 0.05 0.1 0.2 0.5 1.0 




* i i i i i i i i i i i i i i i i i i i i i 
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logiez) 



Figure 5: La relation luminosité-redshift et son application pour mesurer le 
paramètre de décélération q . Les courbes théoriques correspondent respective- 
ment à ço = (en haut), qo = 1/2 (au milieu) et qo = 1 (en bas) obtenues 
en supposant soit A = (courbes continues) soit k = (courbes pointillées). 
Ces mesures sont susceptibles d'être affectées par des erreurs systématiques non 
prises en compte ici (extinction ...) Figure tirée de Perlmutter et al. 1996. 



14 



Facteur 
d'expansion 

redshift 



Changements de phase 
Température 



Echelle de 
Planck 



a ~ exp(Ht) 



a~ t 



1/2 



z=40 000 



a ~ t 



2/3 



z=1100 



t t 



10 A 14Gev 



100 GeV 



1 MeV 



10 eV 



G.U.T., Inflation (?) 
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Découplage des 
photons 



Formation des 
grandes structures 



Figure 6: Histoire thermique de l'Univers. Si les dernières étapes sont relative- 
ment bien comprises, d'autres sont beaucoup plus spéculatives. Par exemple 
il a été envisagé que la baryogenese ait pu avoir lieu aussi bien au moment de 
la brisure de symétrie électrofaible que des énergies de la GUT. Notons que 
la nucléosynthèse primordiale donne une contrainte sur la quantité de baryons, 
fl b h 2 = 0.013 ±0.0003. 
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Figure 7: Le spectre de corps noir observé par l'expérience COBE/FIRAS 
(Mather et al. 1990). 



Me + m p = m + h- ( 50 ) 

Pour les photons on a évidemment /u 7 = 0. Le photon est sa propre antipar- 
ticule et de plus on sait par l'observation du spectre du CMB qu'il n'y a pas eu 
de dépos significatif d'énergie dans le corps noir comme ce pourrait être le cas 
si on avait des particules instables. 

Si on appelle x une espace non-relativiste de masse m x et de dégénérescence 
g x on a, 

fm x T\ 3/2 

n x =g x [-^j-J exp[(/z x - m x )/T\. (51) 

L'équilibre des potentiels chimiques implique alors que (on suppose mn — m p 
dans le pré- facteur) , 



i e n p 1 ( m,,. T 



3/2 

exp[-(m e + m p - mjî)/T]. (52) 



Ce rapport s'exprime facilement en fonction de la fraction de ionisation, x e , avec 

Uenp_ = (53) 
n H n b 1 - x e 
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La température de la recombinaison qui correspond disons à x e — 0.5 a donc 
une dépendance logarithmique avec la densité de baryons parce que le préfacteur 
est très grand: r\ = rib/n 7 w nb/T 3 « 10~ 10 . Sachant que m e + m p — ttlh = 
13.6 eV = 160 000 K, on trouve plus spécifiquement que 

T w 3500X soit z* w 1100. (54) 

Le moment de la recombinaison est donc relativement bien contrôlé. 

A des températures de l'ordre du Mev, ce sont les neutrinos qui se sont 
découplés: il en reste un fond de neutrinos (non-détecté!) dont il faut tenir 
compte quand on calcule la densité de rayonnement. A peu près au même 
moment (à une température légèrement inférieure, tout de même) a eu lieu la 
nucléosynthèse primordiale: formation des noyaux de deutériun, d'hélium 3 et 
4 et de Lithium 7. Dans les conditions de densité de l'univers primordiale, il est 
impossible de former des noyaux plus lourds (ces réactions n'atteignent jamais 
l'équilibre thermodynamique, ce qui donnerait des noyaux autour du fer) parce 
que les noyaux de masse nucléaire 5 et 8 sont instables. La nucléosynthèse 
explosive saute cette étape grâce à une réaction à trois corps, la triple a qui 
donne directement du carbone: cette réaction est évidemment très défavorable 
en milieu ténu. 

Ces étapes sont très bien comprises, et reposent sur de la physique de labo- 
ratoire (les sections efficaces des réactions nucléaires par exemple) mais ce qui 
se passe à température plus élevée est plus spéculatif. Baryogenèse, mécanismes 
de brisure de symétrie sont autant de sujet de recherche actifs. 

Quand on remonte encore en énergie, à l'approche des énergies de GUT 
on a des théories très spéculatives, de plus en plus éloignées de la physique 
de laboratoire. Une idée s'est imposée depuis le début des années 80, c'est le 
mécanisme de l'inflation. 

2.5 Inflation, motivations et principe de base 

L'inflation permet d'expliquer pas mal de chose dans un même scénario. Son 
plus grand mérite est sans doute de mettre dans le champ de la physique une 
problématique générale, les conditions initiales, qui semblait jusqu'à lors en 
dehors de son champ d'investigation. Du même coup l'Univers observable n'est 
plus qu'une infime petite partie de l'Univers au sens plus philosophique du mot. 

• Le problème de la platitude; Aujourd'hui la densité est sinon critique, du 
moins très proche de la densité critique. Cela implique que dans le passé 
la densité était infiniment proche de la densité critique (voir figure 3). 

• Le problème de l'horizon. On observe que des régions du ciel qui, dans une 
cosmologie standard, n'ont jamais été en relations causales ont la même 
température (à 10 -5 près). 



17 




Log(a) 

Figure 8: Représentation schématique du scénario de l'inflation 



• Absence de monopoles. On s'attend naturellement à avoir un reliquat de 
monopoles magnétiques à cause de transitions de phases en théorie des 
Grande Unification. Ce serait des objets massifs qui produiraient une 
densité de 15 ordres de grandeur plus grande que la densité critique. 

• La formation des grandes structures: mécanisme pour la génération des 
fluctuations. 

Comment marche l'inflation? L'idée c'est qu'à un moment donné, après 
l'échelle de Planck, la densité d'énergie est concentrée dans un champ scalaire, 
l'inflaton ip, et que le potentiel de ce champ est à l'origine non nul: cela conduit 
à une croissance rapide du facteur d'expansion. 
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Dans la suite, je ne vais pas essayer de faire une présentation exhaustive de 
la théorie de l'inflation. Je vais m'intéresser à un cas simple correspondant à , 

• un champ simple pour l'inflaton; 

• cas où l'approximation dite de 'Slow Roll' est valable. 
Rappelons que l'action pour un champ quantique s'écrit, 

(55) 




où C est la densité de Lagrangian, 

C = ±{d li <p&><p)-V(<p). (56) 

Le terme supplémentaire qui apparaît dans l'action redonne les équation d'Einstein. 
Dans cette partie, pour reprendre l'habitude de la physique des hautes énergie, 
la constante de gravitation G sera notée, 

G = — L— quand fi = c = 1. (57) 
m P i. 

Dans la suite, je vais illustrer mon propos avec le cas simple, 

V{<p) ex (58) 

correspondant à l'inflation chaotique chère à Linde. Il y a deux dynamiques 
qui sont lices, l'une du facteur d'expansion qui dépend de la densité d'énergie, 
l'autre du champ ip. L'équation d'évolution pour ip s'obtient en cherchant un 
extremum de l'action et on trouve, 

.. à. 1 . dV 

p + 3-p>- Ap = - — . 59 
a a z dp 

L'équation d'évolution du facteur d'expansion s'obtient en écrivant les équations 
d'Einstein. La densité locale d'énergie est donnée par, 

rM = Ç + Œf + vM. (60) 



On a donc, 

87ra 2 \p 2 (Vp) 2 



a 



2 



3«& 



+ V(p) 



k. (61) 



2 2 

En fait d'abord considérer la valeur moyenne du champ qu'on écrit donc, 

p(t,x) =¥>„(*) + (62) 
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et on fait l'hypothèse que la relativité générale s'applique 'classiquement' à 
la valeur moyenne. Cela suppose qu'on s'intéresse à une époque postérieure 
à l'époque de Planck (jusqu'à preuve du contraire on ne sait pas faire de la 
gravité quantique). La quantification du champ se fera sur la partie fluctu- 
ante uniquement, et on suppose que la densité d'énergie portée par <p est bien 
plus importante. Avec cette hypothèse on peut traiter la dynamique de ip 
indépendamment de ses fluctuations. 

L'équation d'évolution donne pour ip Q , 

„ à . d V , . 
<A) + 3-i^o = --3— (ipo), (63) 
a d(p 



et celle pour le facteur d'expansion donne, 

8tt a 2 



à 2 = 



3< 



• 2 

<Po 



V(<po) 



k. (64) 



On reconnaît pour l'équation d'évolution de ip, les termes qu'on attend dans 
une métrique de Minkowski plus un autre terme qui vient de la dépendance en 
temps du facteur d'expansion, qui agit comme un terme de frottement. Pendant 
l'inflation, c'est ce terme qui domine: on a un régime stationnaire où l'énergie 
du champ se dissipe dans une expansion rapide, 

. 9 8tt Go 2 lr , , .„_. 

à 2 = ——V(<p), (65) 

3-ip = -^-(<po). (66) 
a d(p 

La solution de ces équations donnent pour le potentiel de départ, 

<Po(t) = exp (^\f^ m pi -tj > ( 67 ) 

et, 

a(t) ~ exp(Ht), (68) 

avec, 



H(t) = J^f. (69) 



pi. 



Cette solution en a(t) implique que 



2 

a a 



a 



(70) 



ce qui implique une équation d'état effective pour le fluide cosmique p = —p. 
Une conséquence est la rapide décroissance du terme de courbure dans l'équation 
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qui donne la constante de Hubble. En effet flk = —k/a 2 croît pour le big bang 
standard (comme i 1 / 2 ou i 1 / 3 ) mais décroît très rapidement (comme exp(-Ht)) 
pendant l'inflation. 

La validité de l'approximation de 'Slow Roll' implique un certain nombre de 
propriétés de platitute sur la forme du potentiel de l'inflaton. On peut exprimer 
ces conditions comme des conditions sur la forme du potentiel. En particulier 
on doit avoir i/? 2 <C V ce qui implique, 

m pi. fv' x 2 



'-ïSlvl*'. <71) 



et la condition <p'o <C V implique, 



ml V" 



« 1 avec rj = -f^. (72) 

O 7T V 

Cela dit on peut quand même faire marcher l'inflation dans de tels cas. 
Par exemple Starobinsky a examiné un certain nombre de cas où certaines des 
conditions de 'Slow Roll' ne sont pas vérifiées, par exemple quand on a une 
non-analyticité dans la forme du potentiel. 

Mais évidemment l'inflation n'est pas seulement intéressante pour résoudre 
un certain nombre de problèmes lies à l'Univers homogène, cela donne aussi 
un moyen de produire des fluctuations de la métrique à l'origine des grandes 
structures. 

Dans ce cas là, ce sont les fluctuations de l'inflaton que l'on doit examiner. 
L'hypothèse physique sous-jacente est qu'il est légitime de quantifier les fluctu- 
ations de l'inflaton dans une métrique dépendante du temps. Cela implique que 
l'on puisse écrire, 



ôtp 



j d 3 k autpk(t) exp(ik.x) + fâ(t) exp(-ik.x) (73) 



où et a vont être les opérateurs de création de d'annihilation du champ. Ils 
obéissent donc à la relation de commutation suivante, 

[o k ,oi k ,]=<î(k + k / ). (74) 

L'équation du mouvement vérifiée par ce champ est, 

ôifi + 3 H ôip + / =^f = -V"ôip w 0. (75) 

Dans l'approximation du 'Slow-Roll' le terme de source est négligeable et H est 
constant. On a là une équation de champ libre. Cette équation donne la forme 
fonctionelle de ^fe en fonction du temps. 
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Pour avoir la normalisation on peut imposer qu'à petite échelle on retrouve 
les mêmes règles de quantification que pour un espace de Minkowski, c'est à 
dire que 



d r d 



= i<î< 3 >(axi -ax 2 )^ ^ 3 )( Xl - x 2 ) (76) 



Pour une métrique de Minkowski les fonctions ipk sont simplement propor- 
tionclles à des exponentielles et leur amplitude est donnée par 

= Jïjï (^) 1/2 eMikt/a}. (77) 

Pour un espace de Sittcr il faut tenir compte de l'évolution du terme supplémentaire 
en H et on a, 

, / n H ( '. k \ \ ïk 1 

^=(2^( 1+ ^j eXP [^J- (78) 
A petite échelle on retrouve bien sûr les solutions correspondant au cas Minkowski. 
Cela fixe au passage la normalisation de ces modes: c'est parce que les fluctua- 
tions sont à l'origine quantiquc qu'on peut en calculer l'amplitude. A la traversée 
de l'horizon on a un mode croissant très simple. Au bout du compte on a un 
mode simple de la forme, 

» + a -k) . = I d 3 k(^ k e ik x . (79) 

Les conséquences sont nombreuses et très intéressantes: 

• Toutes les observables qu'on peut construire à partir de ces champs vont 
commuter entre elles. On passe donc de fluctuations quantiques à des vari- 
ables qui se comportent comme des grandeurs stochastiques classiques. On 
peut facilement vérifier que les variables aléatoires ^a k + a 1 ^ induisent 
un champ stochastique Gaussien. Pour être plus spécifique on a 

<0\<Pk<Pk'\0> = <5 (3) (k + k')P l p(fc) (80) 
<(# kl ...^ k2p |0> = E il <0|^k,|0> (81) 

combinaisons paires (i,j) 

<0|£ kl ...£ k2p+1 |0> = 0. (82) 

Par la suite on identifie les moyennes sur le vide à des moyennes d'ensemble 
sur les variables 'classiques stochastiques' ip^. 

• C'est la valeur de H qui va fixer les variations de l'amplitude des fluctua- 
tions de densité par la suite. Dans l'approximation du 'Slow roll' la valeur 
de H ne change que très peu pour les échelles accessibles à l'observation, 
et du coup on va avoir un spectre proche d'un spectre Harrison-Zeldovich. 
Pour être plus précis, le spectre de puissance du potentiel est en fc -3 , celui 
de la densité sera en k ns avec n s w 1. 
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• Comme on va le voir par la suite, dans ces approximations on a des fluc- 
tuations dites adiabatiques. 

Conditions pout avoir un modèle d'inflation valide: 

11 faut bien sur que la pente du potentiel soit suffisamment faible pour que ça 
marche. En effet pour que l'inflation remplisse son rôle, i.c. remette causalcmcnt 
en relation des régions de l'Univers observable, il faut que cette expansion ait 
duré suffisamment longtemps. La sortie de l'inflation n'est pas forcément très 
bien comprise: les fluctuations de l'inflaton doivent être prises en compte dans 
la densité d'énergie. On a création de particules. Mais en fait cela n'affecte pas 
forcément le problème de la formation des grandes structures (voir figure [sj) . 

L'amplitude des fluctuations sera de l'ordre de H, donc à peu de chose près 
constante sur les échelles qui nous intéressent: de 1 à 1000 /i _1 Mpc on ne couvre 
qu'une petite fraction du temps d'inflation. 

Cela a un certain nombres de conséquences: 

• Univers à courbure plane; 

• Fluctuations initiales gaussiennes; 

• Un spectre de fluctuation invariant d'échelle dans les cas génériques, c'est 
à dire quand le potentiel de l'inflaton ne présente pas d'accident dans le 
domaine d'énergie qui nous intéresse. 

3 Eléments de cosmologie, l'Univers inhomogène 

3.1 La croissance des fluctuations 

On s'intéresse ici à la croissance des fluctuations dans différents régimes, aussi 
bien pour des échelles plus grandes que l'horizon, plus petites, au moment de 
l'ère du rayonnement que de l'ère de la matière. Les calculs présentés dans cette 
section ont été fortement inspirés de la thèse de Wayne Hu (1995). 

D'une manière générale il faut se donner une métrique qui puisse décrire les 
fluctuations scalaires. On est confronté au problème du choix des jauges, c'est 
à dire qu'il est possible décrire les mêmes fluctuations physiques de différentes 
façons. Pour les fluctuations scalaires on peut choisir la jauge suivante (jauge 
Newtonienne ou longitudinale, voir Bardeen 1980), 

5oo = -a 2 [l + 2V(i)Q(x)] (83) 
gm = (84) 
9ij = a 2 [l + 2 0(t)Q(x)] 7y (85) 

Ici la composante décrit le temps conforme rj avec, 

dt = adt], (86) 
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de telle manière qu'on puisse factoriser a dans la métrique. On va s'intéresser 
à l'évolution d'un certain mode de fluctuation de vecteur d'onde k, c'est à dire 
d'un mode propre du Laplacien à 3D, 

V 2 Q = 7îJ Q = -fc 2 Q. (87) 

Ce mode propre est un champ scalaire. Il définit un champ vectoriel Qi et un 
champ tensoriel Qij définis par 

Qi = Q,i (88) 

Qij = ^Q\v + pijQ- (89) 

L'idée est alors de linéariser les équations d'Einstein en <f> et ip en supposant 
qu'on a des fluctuations de densité dans chacun des fluides cosmiques. 

Evidemment on peut compliquer à l'infini les modèles cosmologiques en in- 
troduisant des tas de composantes, mais pour se fixer les idées je vais sup- 
poser qu'on a essentiellement 3 composantes, des photons (et des neutrinos), 
des baryons et de la matière noire. Cette matière noire sera de plus supposée 
massive, autrement dit ce sont des modèles de type CDM que je vais considérer 
en priorité. 

La partie la plus difficile des calculs à mettre en œuvre concerne la description 
du couplage des photons. On doit écrire l'équation d'évolution de la densité de 
photons /(x, p) dans l'espace des phases. Formellement cela s'écrit sous la forme 
d'une équation de Boltzmann, 

^=cm (90) 

où C[f] est le terme de collisions, diffusion Thomson des photons sur les électrons 
libres. 

Il faut noter que, la diffusion Thomson et les effets Doppler n'affectent pas 
la nature thermique de la distribution d'énergie des photons. Cela justifie que 
l'on identifie les fluctuations d'énergie des photons avec les fluctuations locales 
de température (à un facteur 4 près). On définit, 

-(- 

C'est la température des photons au point x dans la direction 7. L'équation des 
géodésiques s'écrit, 

auquel il faut rajouter un terme de collisions. On peut exprimer le résultat sous 
forme d'une hiérarchie sur les coefficients d'une décomposition multipolaire de 

e, 

0(x, 7 ) = J2 ft(k)P|(k,7)e ik - x (-i)'. (93) 



( X >7) = 7 (Z2- IP 3 àpf-l). (91) 
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L'équation des géodésiques avec terme de collisions donne alors, 



n = 



k 

8 2 = k 
61 = k 



5 



r 1 v 3 



n e x e a T (Oi - V b ) 



— n e x e a T 6 2 



(94) 
(95) 

(96) 



21- 1 



01-1- 



1 



21 + 3 



— n e x e ut Oi pour l > 2. (97) 



On reconnaît un terme de couplage entre les photons et la matière baryonique 
où n e est la densité délectrons, x e la fraction d'ionisation et <tt est la section 
efficace de la diffusion Thomson. 

La partie sans couplage correspond à la diffusion libre. C'est le système 
d'équations qu'on a pour les neutrinos. 

Si je veux pouvoir écrire les équations d'Einstein je dois, pour chaque espèce 
x, écrire le tenseur énergie impulsion, 



T ° = -[1 + <**(*) Q(x)K(t) 
T° = \p x (t)+p x (t)]V x (t)Qi(x) 



T) = Px (t) 



5j + ^3-Q(x)+n x (t)Qj(x) 



(98) 
(99) 

(100) 



On peut alors écrire les équations de continuité pour chacun des fluides jusqu'à 
l'ordre 1 par rapport aux fluctuations de densité. A l'ordre on obtient, 



px = -3 H p x 

V Px) a 

Cette une équation déjà connue. 

A l'ordre 1 des perturbations on a, à partir de T.^° = 0, 



f 



(101) 



Px) \ ' a \p x 



Si on applique cette relation aux photons on a, 
et si on l'applique à la matière, 

<5mat. = - (k Vmat. +3 0^). 



(102) 



(103) 



(104) 
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Pour les photons on retrouve la première équation de la hiérarchie précédente. 
Cela implique en particulier que, 



«mat. - T&y = ~ Knat.)- (105) 



_d_ 

d?7 

On a des équation du mouvement (d'Euler) pour chacun des fluides s'obtiennent 
à partir de T v v l = 0. Pour les photons l'équation est déjà écrite. En appliquant 
cette contrainte au mélange photons+baryons on a 

V b = --V b + k^P + ^n e x e a T (V^-V b ). (106) 

Il reste pour la composante matière froide, 

Vcdm = —VcDM + kip. (107) 

Pour clore mon système d'équation je dois écrire les équations de Poisson 
(i.e. les équations d'Einstein) qui donnent, 

k 2 (V> + 0) = -8 7rGa 2 p T n T ; (108) 

-ij)-4> = AnGa 2 (p T + Z Pt )Y?L ] (109) 
a k 

k 2 4> = ^Ga 2 pr (ôt + s(^) (l+ PT /p T )^j , (110) 

ce qui donne un système fermé. 
Quelques remarques: 

• Les indices T font référence à des quantités totales (obtenues en sommant 
sur tous les fluides pondérés par leur densité moyenne respective). 

• Ces équations dépendent de la jauge utilisée. En particulier le terme de 
source de l'équation de Poisson proprement dite (la dernière) s'identifie 
dans une certaine jauge à la fluctuation de densité totale. 

• IIt est la partie anisotrope de la pression. Elle est nulle pour la matière 
(p = de toute façon). Elle est nulle pour les photons quand ils sont 
couplés aux baryons. La seule contribution résiduelle vient des neutrinos 
après qu'ils se soient découplés de la matière. 

A partir de là je peux chercher à calculer les différents modes de croissance. 
Le schéma de la Fig. 9 montre les différents régimes et leur enchaînement. 

Une quantité importante est la traversée de l'horizon. Dans l'espace de 
Fourrier cela est donné par kî] = 1 (77 est le temps conforme). 
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COBE/DMR 1111111 m H Cat. de galaxies 

Planck WMMWMMmMMm 

Figure 9: Les différents régimes de croissance gravitationnelle. A donne la 
croissance du contraste de densité de la matière noire froide. Les zones hachurées 
du bas montrent les différents domaines accessibles aux observations. 
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Figure 10: Evolution de l'amplitude du potentiel gravitationnel à grande échelle 
au moment du passage de l'équivalence. La courbe du haut donne la variation 
de é (source de l'effet ISW) et la courbe du bas, 6 + 4>. 

• A une échelle plus grande que l'horizon, krj <C 1, Univers dominé par le 
rayonnement. On a 



d r r \ à 

drj 



(Vmat.-Vg = — 0W -F 7 ); (111) 



^(<W-^ 7 ) = fc(v mat .-y 7 ). (ii2) 

Dans une théorie inflationnaire à champ simple la matière tombe dans les 
potentiels ensemble. Il n'y a pas de création d'entropie, donc, 

3 

£>mat. = |&y, (H3) 
U m at. = Vy. (114) 

Dans le cas contraire (inflation hybride...) on peut avoir des fluctuations 
isocourbes: fluctuation de la composition de l'Univers sans fluctuation de 
courbure. 

Finalement on doit résoudre, 

^ + ^ = -\(^) 5 ^ ^7-4^, (115) 

ce qui donne après avoir éliminer ô 7 

d 2 tb dé 

-4+4-f=0 116 
der da 
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dont la solution générale est, 

ip~tpc + i>D\- (H7) 
cr 

On a un terme 'croissant' et un terme 'décroissant'. Pour ce terme dit 
croissant. Les fluctuations de densités (que j'identifie de manière générale 
avec le terme de source de l'équation de Poisson) croissent comme k 2 ï] 2 . 

La traversée de l'équivalence à une échelle plus grande que l'horizon. 
L'équation générale s'écrit en fait, 

à , ; 1 f à\ <L 4 + 3a/a eQ i , ; ,„„^ 
a Y 2 V»/ 4 l + a/a cq . K ' 

Il existe une solution explicite (voir Fig. ). Finalement on a, 

0) (119) 
oo), (120) 

et on remarquera que <5 7 /4 s'identifie à la fluctuation locale de température. 

La traversée de l'horizon. C'est là que les choses deviennent plus com- 
pliquées puisqu'il faut tenir compte des intéractions entre photons et 
baryons. Tant qu'on est avant la recombinaison ces 2 fluides restent 
couplés, mais la matière non-baryonique n'est plus couplée à ceux-ci. 
Quand le rayonnement domine on a, 

<5 7 + y<5 7 = 4^-^ 2 V. (121) 

Du coup les fluctuations sont piégées dans des ondes acoustiques. La 
matière non-baryonique est spectatrice. Quand le rayonnement domine 
l'équation d'évolution des fluctuations de matière sont alors, 

6cDM + -6 C DM = -&î>-k?il;-3-4>. (122) 
a a 

Le membre de droite est un terme de source déterminé par le comporte- 
ment des photons. Comme les croissances des fluctuations pour les pho- 
tons sont interrompues, on a essentiellement une équation homogène dont 
la solution croissante est une croissance logarithmique de la fluctuation de 
densité. 

Univers dominé par la matière. Cette fois le terme de source de l'équation 
d'évolution est donné par la matière elle-même. En fait je retombe sur 



ip(a — > oo) 
<5 7 (a — > oo) 



10 



V>(a 
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une équation déjà écrite précédemment, et les potentiels sont de nouveau 
constants: les fluctuations de densité croissent comme rf <~ a. 

Par contre le fluide photons+baryons connaît toujours des oscillations 
acoustiques. La vitesse du son c s est maintenant donnée par, 

Cl = 3(l + p b /p 7 )- (123) 



On a essentiellement, 



5, 



7 



cos(fc / c s d?7) (124) 



-<k r 

Jo 

r 

sin(fc / 
Jo 



V 7 ~ sin(fc / c s dr/). (125) 
Jo 

(si c s était constant cela ferait kr\j\fï dans l'argument des fonctions 
trigonométriques) . 

Quelles sont les conséquences observationnelles? Il y a essentiellement 2 as- 
pects qui sont importants d'un point de vue obscrvationncl. Les fluctuations de 
températures sur le 3K et les fluctuations de densité qui vont données naissance 
aux grandes structures de l'Univers. 

• Les fluctuations de température. Localement on peut décomposer les fluc- 
tuations de température du rayonnement dans son espace des phases par, 



9{j) s / JVdp/(x, P ) _ A /4 = E aim Yr{il (126) 

où 0(7) est la fluctuation de température dans la direction 7, décomposée 
en harmonique sphérique (implicitement à la position de l'observateur). 
Pour un champ de fluctuations isotropes on a 

< m <H>,m' >= àlV ^mm' C; ■ (127) 

Cela définit les Ci. 

Les anisotropies du 3K permettent de voir en quelque sorte l'état du 
plasma photons+baryons au moment de la recombinaison. Après la re- 
combinaison les photons se découplent complètement, leur libre parcours 
moyen devient infini, ils diffusent donc librement. Les anisotropies de 
températures sont alors la superposition de 3 effets différents: les fluctu- 
ations de température intrinsèques, les fluctuations de potentiel qui in- 
duisent des effets de redshift gravitationnel, enfin des effets Doppler dûs 
aux mouvements du fluide le long de la ligne de visée. Il y a un cut-off 
naturel aux petites échelles dû au fait que la surface de dernière diffusion 
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Figure 11: Evolution de l'amplitude des fluctuations intrinsèques de température 
(en haut) et de la vitesse particulière (en bas) en fonction du facteur d'expansion 
et du mode k. 
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Figure 12: Exemples de spectre de fluctuations de température pour des modèles 
standards CDM. C'est la quantité Ci 1(1 + 1) qui est donnée en fonction de l. 
Différentes valeurs de ÇIq h 2 sont proposées (figure tirée de Hu, 1995). 



a une épaisseur finie (ce qui induit un 'Silk damping'). L'amplitude du 
mode l est donnée essentiellement par, 



Ôi 



21 + 1 



(OoM + iPMMkv*) + 9 ^l^( k v)+ (128) 



I/O 



exp(-fe 2 /fc2), (129) 



pour un mode k. Cette expression est valable dans la limite d'une surface 
de dernière diffusion infiniment mince. Le première correction à cette 
limite donne un effet d'amortissement à petite échelle, le 'Silk damping' à 
une échelle ko liée à la longueur de diffusion des photons. En général il 
faut résoudre un système d'équations couplées pour obtenir les différents 
moments. La projection des fluctuations locales de température sur les 
harmoniques sphériques fait intervenir des fonctions de Bessel ji pour le 
cas d'un univers plat. Pour un univers avec k ^ 0, c'est plus compliqué. 
On a alors un spectre C; typiquement de la forme donnée par la figure ^2|. 

A ce jour, on ne connaît correctement que la partie à grande échelle grâce 
aux observations de COBE/DMR. La forme des pics acoustiques est une 
chose que les gens cherchent activement à mesurer parce qu'ils contiennent 
des informations précieuses sur les paramètres cosmologiques. Il y a en 
préparation de nombreuses expériences sol ou ballon, sans compter les 



32 



projets satellitaires américains avec MAP ou européens avec Planck, qui 
doivent mesurer ces anisotropies avec des sensibilités et des précisions 
diverses. 

• La conséquence observationnelle qui va davantage m'interesser ici est le 
spectre de puissance P(k) des fluctuations de densité de la matière. Dans 
les modèles avec de la matière noire en quantité dominante les baryons ne 
vont contribuer à la forme de ce spectre que de manière mineure. 

Du coup fc C q , échelle qui traverse l'horizon au moment de l'équivalence 
qui est un paramètre essentiel de la forme du P{k). En effet les modes 
k inférieurs à fc eq . croissent régulièrement, comme t] 2 , alors que les modes 
avec k plus grand que fc cq . sont gelés entre la traversée de l'horizon et 
l'équivalence. 

Pour les modèles les plus réalistes, le P(k) a une forme simple sans struc- 
tures particulières. En particulier les oscillations plasma ne se voient que 
lorsque la fraction de baryons est suffisante. 

En résumé on a 

P{k) ex k pour k > fc eq ., (130) 
P(k) ex k~ 3 pour fc<fc cq ., (131) 

pour un spectre initial Harrison-Zeldovich. 

3.2 Pourquoi un modèle avec de la Matière Noire Froide 

? 

La première raison est de rendre compatible le niveau des anisotropies du 3K 
avec les fluctuations de densité dans l'Univers local. Au moment de la recom- 
binaison les fluctuations de densité de la matière noire peuvent en effet être 
beaucoup plus grandes que celles des photons+baryons. Les baryons tombent 
dans les puits de potentiel de la matière noire après la recombinaison. Du coup, 
on peut très bien avoir des anisotropies de température à un niveau bas tout en 
ayant des fluctuations de densité de la matière noire relativement importantes. 

3.3 Évidences de l'existence de matière noire 

Evidemment cette belle idée doit être soutenue de manière expérimentale en 
essayant de détecter la matière noire non-baryonique de manière directe. 
Les évidences sont les suivantes, 

• Les courbes de rotation des galaxies. Cela s'applique aux halos des galax- 
ies. 
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Figure 13: Evolution à petite échelle du contraste en densité des baryons 
comparé à celui de la matière noire froide. Les baryons restent couplés très 
longtemps au photons. Au moment de la surface de dernière diffusion, le con- 
traste en densité atteint par la matière noire froide peut être beaucoup plus 
important que celui des baryons (figure tirée de Hu, 1995). 



• Les études dynamique des amas de galaxies. Un amas de galaxies est un 
immense puit de potentiel dans lequel les galaxies sont liées. La mesure 
de leur dispersion de vitesse oermet d'estimer la masse de l'amas. 

• Les mesures de masse par effet de lentille. Exemple A2218 avec un système 
d'arcs. C'est très propre puisque c'est une détection directe de la masse 
projetée, par d'hypothèse sur l'équilibre hydrodynamique dans l'amas. 
C'est une méthode très prometeuse car elle peut s'appliquer à toutes les 
échelle, aussi bien pour sonder les halos des galaxies, que pour faire des 
détection de la matière gravitante à grande échelle. 

• Les champs de vitesse à grande échelle. L'idée est de mesurer les vitesses 
particulières, écarts des vitesses d'éloignement avec le flot de Hubble, à 
l'aide d'indicateurs spécifiques (de type dispersion de vitesse). C'est un 
domaine potentiellement très intéressant, cependant il souffre d'un défaut 
important: les erreurs de mesure augmentent avec la distance. Au delà de 
50 h~ 1 Mpc les erreurs dépassent les 100% sur une galaxie individuelle. 
Cela vient des dispersions intrinsèques des relations utilisées et pas des 
erreurs de mesures. 
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3.4 Modèles alternatifs 



Les modèles avec défauts topologiques sont de moins en moins crédibles en tant 
que mécanisme pour former les grandes structure. La principale raison est qu'il 
est très difficile de concilier le niveau des anisotropies vues par COBE avec 
l'amplitude des fluctuations à grande échelle (voir l'article récent de Albrecht 
1998). 

Le calcul de la croissance des fluctuations dans de tels modèles est rendu 
très difficile par l'existence de termes de sources dans les équations d'Einstein. 
Il faut tenir compte de la dynamique propre des sources, et des corrélations 
non-triviales qu'elles induisent à des temps différents. 



4 La dynamique gravitationnelle, les théories 
linéaires 

A partir de maintenant je m'intéresse à la dynamique de la matière noire en 
effondrement gravitationnel. 

4.1 Equation d'évolution dans l'espace des phases 

Sous l'horizon l'équation dynamique s'écrit, 

du à f ,» ,J(x')(x'-x) ,.,„„. 

- + -u = ë = G p flj /dV yi_ x|3 > , (132) 

où u est la vitesse particulière des particules (écart de la vitesse totale avec le 
flot de Hubble). L'impulsion s'identifie à (c'est aussi la variable conjuguée de x 
dans le Lagrangien décrivant la dynamique), 

p = uma. (133) 

J'écris la densité dans l'espace des phases, 

/(x,p)d 3 xd 3 p, (134) 

où x est la position comouvante et p est l'impulsion. Le théorème de Liouvillc 
appliqué à l'évolution de la densité / donne, 

f = m f <* p ' t] + iè /(x ' p ' t] - mV -* (x) l /(x ' p ' t] = (135) 

A$(x) = ^|/(x, P ,t)d 3 P . (136) 
Dans toute sa généralité cette équation est évidemment très difficile à résoudre. 
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On peut prendre différents moments de cette équation par rapport à p. La 
densité de matière est par construction, 

p(x,i) = ^ J d 3 p/(x,p). (137) 

J'écris cette densité 

p(x,t)=p(t)[l + 6(x,t)], (138) 
et alors p(t) oc a 3 . La vitesse moyenne est u donnée par, 

u= Jd'pp/(xp) 
ma J j d 4 p 

et je vais écrire de manière un peu cavalière que, 

Jd 3 pp iPj /(x,p) 

m 2 a 2 //d3p u^ + ^u^. il4Uj 

Quand on intègre l'équation de Liouville par rapport à p on obtient, 

pf t + ^V.p(x) u(x) = 0, soit , ^ + ±V.(1 + <5(x))u(x) - 0, (141) 

c'est à dire une équation de continuité. 

Quand on intègre l'équation par rapport à p après l'avoir multiplier par p 
on a, 

^ J Pi fd 3 p+-^dj ; Jp iPj fd 3 p + a 3 p(x,t)$ i = 0. (142) 
Je peux réécrire cette équation en faisant apparaître la vitesse particulière, 

W + h + = 4 Vi$ - Ta • ( W3 ) 

Faisons des hypothèses supplémentaires. Si localement le fluide est ther- 
malisé on a 

{muj) c = 6i^. (144) 

La dispersion de vitesse locale s'identifie avec la pression en supposant qu'elle 
est isotrope. On a alors, 

dui à 1 . „ . 1 , Pi ,„ ,„n 

_l + -u i + - (uiVj )uj =~Vj0-^. 145 
dt a a a pa 

Si on se donne une équation d'état on peut alors fermer le système. 
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4.2 La longueur de Jeans 

Une question qu'on peut se poser est de savoir à quelle échelle la pression peut 
jouer un rôle et freiner les effondrements gravitationnels. Pour faire ce calcul il 
faut se donner une expression pour la pression. 
Par exemple, pour un gaz mono-atomique on a, 

£ = ^ (146) 
dp 3 m 

On peut ensuite linéariser les équations d'évolution, ce qui donne 

^ + -V.u = (147) 
ut a 

puisque 

p = p + c 2 s pô. (149) 
Cela définit une longueur caractéristique, la longueur de Jeans, avec 

A, 7 =c s (VGp) 1/2 . (150) 

Les fluctuations de taille plus grande que Aj s'effondrent librement sans être 
contrer par la pression. Cela donne des échelles de masse de 10 5 M©. Pour les 
échelles cosmologiques (i.e. la masse d'une galaxie est typiquement de 10 12 M©, 
celle d'un amas de fO 15 M©) on peut donc en toute sécurité négliger la pression 
due aux fluctuations thermiques. 

Notons quand même qu'admettre que la pression est nulle c'est dire que 

(u iUj ) c = (151) 

donc la densité dans l'espace des phases peut s'écrire 

/(x,p) =/)(x)i (3) [p-mdu(x)]. (152) 

Cela revient à supposer qu'en tout point on a un seul flot. Cela ne sera vrai qu'au 
début de la dynamique. Ultérieurement on s'attend à avoir des croisements de 
coquilles (voir fig |Î4|) . 

4.3 L'approximation Newtonnienne avec un seul flot 

Dans l'approximation à un seul flot, on a finalement le système d'équation suiv- 
ant, 
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Figure 14: Description schématique de l'espace des phases après les premiers 
croisements de coquilles. 



^p(x,t) + -V x .[p(x,i)u(x,i)] = 
ot a 

Q ' 1 

— u(x,t) + -u(x,t) + -(u(x,t).V x )u(x,t) = 
Ot a a < 

V x 2 ^(yL,t)-4nG[p(x,t)-p(t)]a 2 = 0. 



(153) 

--VxV(x,i) (154) 



(155) 



C'est ce qu'on appelle une description Eulcriennc de la dynamique. 



4.4 Description Eulérienne ou description Lagrangienne 

On peut en faire une description équivalente en ayant une approche Lagrangi- 
enne. L'idée est de décrire la trajectoire des particules à partir de leur position 
dans l'espace de départ q. On écrit donc qu'une particule initialement en q 
se trouve à un instant t à la position (comouvante) x après s'être déplacée de 
*(q,t), 

x = q+#(q,i). (156) 

Le champ déplacement ^(q, t) est donné par la même équation d'Euler sim- 
plement la densité est donnée par l'inverse du Jacobien de la transformation de 
q vers x. Plus explicitement on peut montrer que 



<9x 1 



(V q .*) 2 -£*iJ*i,i 



(157) 
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+è [( v q -*) 3 - 3V q -*E« + 2E ijfe *^,**m 

On va voir que selon le problème qui nous intéresse on pourra prendre l'une 
ou l'autre des formulations. 

4.5 La croissance des fluctuations en théorie linéaire 

La théorie linéaire consiste à supposer que les fluctuations sont de faible ampli- 
tude, ô <C 1. Pour une description Lagrangienne, cela veut dire que les gradients 
du champ de déplacement sont faibles, et donc que le Jacobien de la transformée 
est proche de 1. 

Dans le régime linéaire il est plus simple de regarder la description Eulériennc. 
La linéarisation des équations donne, 

| + iv.u = 0, (158) 

jg+2-j£ = (159) 
ot z aot 

La deuxième équation permet de rechercher les modes de fluctuation. On 
trouve 2 modes, 

6{x,t) = D + (t)5+(y)+D-(t)5-{x), (160) 

dont la dépendance en t dépend des paramètres cosmologiques (elle est cachée 
dans la dépendance en temps du facteur d'expansion). Pour être plus spécifique 
on a, 

D + 2HD- ^H 2 (t)fl(t)D = (161) 

qui est valable pour tous les modèles cosmologiques. 
Pour un Univers Einstein-de Sitter on a, 

D+(t) ex i 2/3 , (162) 
D_(t) ex l/t. (163) 

Notons que D + peut s'exprimer analytiquement en fonction de t si A = 0. Il 
doit être calculer numérique pour des modèles avec constante cosmologique. 

Une quantité utile est la dérivée logarithmique de D + avec le facteur d'expansion, 
puisqu'elle apparaît de manière naturelle dans l'équation de continuité. On a 
un fit analytique de la forme (Lahav et al. 1991), 

m A) = = n °' e + tïï a (1 + n/2) - (164) 

On voit que ça dépend très peu de la constante cosmologique. 
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Une relation intéressante en théorie linéaire est la relation densité-vitesse. 
Localement elle s'écrit à partir de la relation de continuité, 

0(x) = J-V x .u(x) = -f(il, A)5(x). (165) 
a II 

Remarquons que est bien a priori une quantité observable qui ne dépend pas 
de la constante de Hubble. Les flots à grande échelle sont connus en km/ s, 
à des distance qui sont déterminées en km/ s. Plus précisément 9 exprime les 
fluctuations locales de la constante de Hubble, i.c. du taux d'expansion. 

Cette relation a évidemment une contrepartie non-locale qui exprime la re- 
lation entre la vitesse locale et les fluctuations de densité environnantes, 

JL = m [ M{ ^^Lz2L_ (166) 

a H 4ir J |x-x'| 3 

Cette relation exprime que la vitesse acquise par une particule en régime linéaire 
est proportionnelle à son accélération. 

A ce niveau on peut facilement donner une interprétation des modes crois- 
sants et décroissants. Le mode croissant correspond au cas où le champ de 
densité et de vitesse sont en 'phase': quand la matière tombe vers les puits 
de potentiel. On peut très bien imaginer qu'on puisse construire un champ de 
vitesse qui soit tel que que les particules aient tendance à s'échapper des puits 
de potentiel et tendent ainsi à effacer les fluctuations initiales. Ce serait le mode 
décroissant. 



4.6 L'approximation de Zel'dovich 

Cela consiste à prendre la solution linéaire pour le champ de déplacement mais à 
garder l'équation de continuité dans toute sa généralité (Zel'dovich 1970). Dans 
une telle approximation, le champ de déplacement s'écrit comme la somme de 
deux termes, 

^(q,t)=£> + (t)^+(q)+D_(t)V-(q), (167) 

où D + (t) et D-(t) sont les fonctions du temps solution de l'équation précédente. 
On retrouve des termes croissants et des termes décroissants. Quand on ne garde 
que le terme croissant, on a une factorisation dans la dépendance en temps et 
la dépendance en espace. Les trajectoires des particules sont donc rectilignes et 
suivent la direction de la force appliquée aux particules à l'instant initial. 

La densité locale peut alors s'exprimer en fonction des valeurs propres de la 
matrice de déformation. Plus précisément on a, 



(1 - Ai £>+)(l-A 2 1>+)(1- As !>+)' 
où les Ai sont les valeurs propres de la matrice tpij. 
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Pour des conditions initiales Gaussiennes on peut alors calculer la distri- 
bution de densité locale à partir de la distribution des valeurs propres. La 
distribution des valeurs propres est donnée par (Doroshkevich 1970), 



p(Ai, A 2 ,A 3 )dAi dA 2 dA 3 = dAi dA 2 dA 3 



5 5 / 2 27 

87TCrg 



(Ai - A 2 )(Ai - A 3 )(A 2 - A 3 ) exp 



O ( 3Jq^ „ J{ 



0"n 



02 



(169) 
(170) 



ou 



^oi — Ai + A 2 + A 3 , 

J02 = AiA 2 + Ai A 3 + A 2 A 3 . 



(171) 
(172) 



On peut alors calculer la fonction de distribution de probabilité de la densité 
locale dans cette approximation (Kofman et al. 1994), 



P{p)A P = 



9 5 3 / 2 dp 
47r7V s( o 3 cr 4 y 3( l)l/3 



ds e 



-(s-3) 2 /2<t 2 



-m/2, 



Pn{s) = S 5 



1/2 



+ cos 



2 , , 1 

-An — l)7r + - arecos 

O ô 



54 
ps 3 



1 



(173) 
(174) 

(175) 



où paramètre a(t) = D(t)ao est la déviation standard des fluctuations de densité 
p dans la théorie linéaire et N s est le nombre moyen de flots (N s — 1 dans le 
régime à 1 flot). 



4.7 La vorticité 

Il n'est pas inutile de remarquer que la vorticité va être diluée par l'expansion. 
Si on définit la vorticité par, 

w k = (uij-u jti )^ k , (176) 
(e*- 7 '* est le tenseur totalement antisymétrique) on a (en régime linéaire), 

w = — w. (177) 
a 

Il n'y a évidemment pas de terme de source. La vorticité est un mode décroissant: 
elle est diluée par l'expansion (i.e. conservation du moment angulaire). 

Cette propriété est en fait vraie à n'importe quel ordre de la théorie des 
perturbations. (Notons que sa mise œuvre pour l'approche Lagrangienne est 
loin d'être triviale). 
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4.8 Domaine de validité de l'approximation linéaire 

C'est valable au delà de 10 ft, _1 Mpc. En effet la r.m.s. des fluctuations atteint 1 
vers 8/i _1 Mpc. Les expériences numériques montre qu'une description linéaire 
est correcte au delà de cette échelle. On verra par la suite qu'à ces échelles 
on peut faire des choses plus sophistiquées que la théorie linéaire, mais d'une 
manière générale l'approche perturbative est valable à ces échelles. 

Il est important d'avoir à l'esprit que la théorie linéaire ne change pas la 
forme du spectre de fluctuation mais uniquement son amplitude. Pour résumer 
on a, 

P(fc) iocai = 2 p(fc) rccomb - = (§ ± |^y) 2 T 2 (k) P(fc) prim ', (178) 

où T(k) décrit la fonction de transfert qui contient toute la micro-physique de 
la recombinaison. 



4.9 Les galaxies comme traceurs du champ de densité 

Un moyen économique (mais relativement dangereux) de se faire une idée de 
la forme du spectre est de considérer la distribution des galaxies dans notre 
Univers local. 

On a maintenant des données assez précises sur la forme du spectre dans 
notre Univers local. Une forme phénoménologique est donnée par (spectre 
obtenu à partir du catalogue angulaire APM par Baugh et Gaztanaga 1996), 

m - + ("») 

avec 

kc = h h Mvc ~ x - ( 18 °) 

La normalisation du spectre est telle que la variance des contrastes en densité 
dans une sphère de rayon 8/i _1 Mpc est de l'ordre de l'unité, 

tT » = (\T^Ô5 / <5(x)d 3 x) ) 1/2 wl pour R 8 = 8 h' 1 Mpc. (181) 
\4tt.R 8 J| X | <JÎ8 J 

La valeur de as pour les galaxies varie un peu avec la population d'objets 
sélectionnés. On observe en particulier une légère augmentation de ag avec 
la luminosité des objets. 

Pour ces raisons historiques, la normalisation du spectre de P(k) est souvent 
discutée en terme de as- D'un point de vue dynamique il faut retenir que c'est 
une échelle de transition. 
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5 Vers la dynamique non-linéaire: l'effondrement 
sphérique 

L'émergence de non-linéarités est en soit un problème difficile. Un cas simple 
pour lequel on peut calculer l'évolution de manière complète est de regarder 
l'effondrement d'une région à symétrie sphérique. 

Je suppose donc que j'ai une fluctuation initiale de forme, 

^(x) = /(|x-x |). (182) 

Pour simplifier on peut penser à un profil en marche d'escalier mais ce n'est pas 
nécessaire. A cette fluctuation on peut appliquer le théorème de Birkhoff, donc 
on peut écrire l'équation d'évolution du rayon de cette perturbation, 

M « R > = ('-47ë)f R '- < 184 > 

Tant qu'on a pas de shell-crossing, la masse (hors À) se conserve. On a finale- 
ment une équation d'évolution tout à fait similaire à celle du facteur d'expansion, 
mais avec une constante d'intégration différente. La solution de cette équation 
est connue (quand À = 0). Il reste une difficulté technique qui est d'exprimer la 
solution R en fonction de a. Quand = 1, c'est relativement simple (a et t sont 
simplement reliés) mais sinon la relation est formellement assez compliquée. 

Les conditions initiales sont déterminées par la valeur de la sur-densité ini- 
tiale, 

*-ds5î- 1 - (185) 

où po(t) est la densité moyenne de l'univers, et par la taille initiale Ri de la 
fluctuation. 

Dans le cas d'un univers plat sans constante cosmologique, la sur-densité 
4ttMpo ne dépend que du produit de la densité initiale et du facteur d'expansion 
dont l'univers a crû depuis le début de l'effondrement. La sur-densité s'écrit, 
quand la sur-densité initiale est positive, 



3M 9(0-sinl 



4irR(t) 3 p (t) 2(1-008 0)3 

et, quand la sur-densité initiale est négative, 

o(t) 3 /3 



(187) 



a(U) 5 V4 



2/3 

(sinh 6-6)) (188) 
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3M _ 9(sinhtf-0) 2 



4irR(t) 3 Po (t) 2(cosh6»-l) 3, 

Il est intéressant de remarquer que la densité (exprimée comme le rapport 
M < (R)/R 3 ) est donnée en fonction de la fluctuation linéaire Si D + approxi- 
mativement par, 



Cette relation (exacte quand Çl — ► 0) est valable pratiquement pour n'importe 
quelles valeurs des paramètres cosmologiques. Ce qu'on voit c'est qu'au bout 
d'un temps fini, on rencontre une singularité, i.e. la densité devient infinie. Cette 
singularité est atteinte quand Si D+(t) — 5 C avec ô c — 1.69 pour la dynamique 
exacte dans un espace Einstein-de Sitter (dans la dynamique approchée ci-dessus 
ô c = 1.5). Si Si est négatif on ne rencontre pas de singularité, et la densité se 
comporte comme, 

p~ {-ôiD + )-V 2 , Si<0. (191) 
Les régions initialement sous-denses se vident donc peu à peu. 

5.1 La virialisation 

La description de l'effondrement par les équations précédentes n'est valable que 
lorsqu'on n'a pas de croisement de coquilles. Dès que l'on a des croisements, la 
masse M(< R) n'est plus conservée. 

Ce qui va se passer alors est un processus de virialisation, la matière va 
tendre vers un équilibre dynamique. Pour évaluer l'état final du système on 
peut faire un simple bilan énergétique. Pour un système isolé on a l'équation 
suivante, 

avec, 

K ' S (193) 

particules 

y - E- o Gm2 r ( 194 ) 

^ 2ax,-x, y ' 

paires J 

Si on suppose que l'équilibre est atteint quand le second membre s'annule, alors 
le rayon de la structure qui se forme ainsi doit être la moitié du rayon de 'turn- 
around', rayon physique (et non-comouvant) maximal atteint par la structure. 

Cette idée a été vérifiée dans des simulations numériques. Au moins pour 
des structures suffisamment isolées, c'est une bonne description de ce qui se 
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passe. Le contraste en densité qui en résulte A, est de l'ordre de 200, 

200. (195) 

Les mécanismes de virialisation sont cependant loin d'avoir été élucidés. 



^ _ Pobjct. 

Po 



5.2 La théorie de Press et Schechter 

Comme application du modèle du collapse sphérique, on peut chercher à calculer 
le nombre d'objets virialisés qui sont formés à une époque donnée. L'idée est la 
suivante, on cherche la fraction de matière qui se trouve dans un objet virialisé 
de masse donnée. 

La masse d'un objet par conservation de la masse est donnée par la taille R 
avec M = 4ir R 3 /3 de la perturbation initiale. On va donc chercher la probabilité 
qu'un point donné de l'Univers se trouve initialement dans une fluctuation de 
densité à l'échelle R qui soit telle qu'elle est virialisée à une époque ultérieure 
donnée. Cette probabilité donne alors la fraction de masse de l'Univers /(> M) 
qui est dans des objets de masse plus grande que M, 



f(>M)= f àSiP{Si,R) 

J5c D + (ti)/D+(t ) 



(196) 



où p(ô) est la fonction de distribution de probabilité de la sur-densité initiale 
Si à une échelle de filtrage R. Pour des conditions initiales gaussiennes p est 
simplement une gaussienne de variance Oi(M) qui s'exprime en fonction du 
spectre et de la fenêtre de filtrage. Alors, 



/(> M) = \ ( 1 - Erf 



Ô C D + (U 



V2D+(t )ai(M) 



(197) 



et donc 



/(M)dM 



-S c 



D+(to) 



<Ti(M) 



dlogMM)) 
dM 



exp 



D+(to) 
D+(U) 



af(M) 



dM. 



D+(ur , 

" (198) 

C'est une relation séduisante mais par bien des égards, très suspecte. On peut 
déjà remarquer que la fraction totale de matière dans des objets virialisés atteint 
1/2 alors qu'on s'attendrait à ce que cela atteigne 100%. Il se trouve qu'en 
corrigeant cette formule d'un simple facteur 2 on a un bon accord avec les 
résulats de simulations numériques. 

Le problème est que ce calcul ne prend pas en compte toutes les structures 
qui sont finalement susceptibles d'entraîner un point de matière donné. 

Bond et al. 1991 ont donné une interprétation alternative à cette formule en 
prenant comme fenêtre de filtrage un filtre carré dans l'espace des phases. Alors 
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à chaque changement d'échelle la densité locale est donnée par une quantité 
supplémentaire qui est indépendante de la densité locale à l'échelle précédente. 
Le problème finalement peut se formuler sous forme d'une marche aléatoire. 
Cette approche a le mérite de redonner le facteur 2 (toute marche aléatoire 
croise une frontière fixe à distance finie), mais il reste que le problème n'a pas 
été résolu dans toute sa complexité (voir Blanchard et al. 1992). 

En plus de la complexité du problème purement statistique, la dynamique 
d'un élément de volume pris au hasard n'est certainement pas toujours locale, 
i.e. dépendante uniquement de la surdensité locale, indépendante des effets de 
marées. 

La théorie de Press et Schechter est en tout état de cause plus pertinente 
pour des objets isolés, rares, comme les amas de galaxie. 

5.3 La densité d'amas pour contraindre l'amplitude du 
spectre de puissance P(k) 

La densité des amas de galaxies est très utile pour contraintre l'amplitude du 
spectre de fluctuation, puisqu'elle teste le comportement de la distribution de 
densité initiale dans sa coupure exponentielle. Typiquement un amas correspond 
à un événement à 2, 3 a par rapport aux fluctuations initiales. Même si il y a 
quelques incertitudes sur les préfacteurs, cela n'affecte pas trop les conclusions 
qu'on peut tirer sur la valeur de l'amplitude des fluctuations donc de as ■ 

Comme on voit que la quantité qui est effectivement testée est une combinai- 
son de D+ et de a on voit bien que l'on va contraindre une combinaison de as 
et des paramètres cosmologiques (avec une dépendance sur la forme du spectre 
que je vais ignorer dans la suite). 

En pratique on a (voir par exemple Oukbir & Blanchard f997, Eke et al. 
1996), 

as^o 5 ~ 0-6 ±0.15. (199) 

Cette analyse dépend un peu de la forme du spectre. L'échelle de masse des 
amas est bien 8 ft, _1 Mpc pour un Univers Einstein-de Sitter mais ce n'est pas 
forcément vrai autrement. Notons que la dégénérescence entre Clo et as peut 
être levée en examinant la dépendance en z de la distribution des amas (voir 
Oukbir et al. 1997). 

6 Le régime quasilinéaire, effets des couplages 
de modes 

Dans cette partie je vais explorer la théorie des perturbations appliquée à la 
formation des grandes structures, d'un point de vue théorique mais aussi en 
illustrant ces techniques par des applications observationnelles possibles. 
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Figure 15: Evolution quasi-linéaire de la forme de la fonction de distribu- 
tion de probabilité P quand on entre dans le régime non-linéaire. L'évolution 
différentielle entre les parties sur-denses et les parties sous-denses induit une 
asymétrie. 

Un des objectifs de ces calculs est de déterminer des quantités statistiques 
qui tiennent compte aussi bien des propriétés statistiques initiales que de la dy- 
namique. Une quantité plus particulièrement considérée (parce qu'elle est facile 
d'accès observationnellement aussi bien que numériquement) est la fonction de 
probabilité de la densité locale, P(S). 

En pratique le champ de densité est obtenu après filtrage d'une représentation 
discrète de ce champ, par des points d'une simulation numérique ou des galax- 
ies dans un catalogue. Alors P(S)d5 est la probabilité que la densité locale soit 
entre ô et dô. 

Initialement on sait que cette distribution doit être gaussienne. La largeur de 
cette distribution dépend de l'amplitude et de la forme du spectre de puissance 
P(k). Plus précisément, 

° 2 {R) = J j^P(k)W 2 (kR), (200) 

où R est l'échelle de filtrage et W est la fonction de filtrage dans l'espace des k. 

Les résultats obtenus sur l'effondrement sphérique donnent une idée de ce qui 
peut se passer quand les non-linéarités commencent à jouer un rôle: l'effondrement 
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des régions surdenses est accélérés alors que les régions sous-denses se vident de 
moins en moins vite. La figure ( |Î5| ) donne une idée de ce qui en résulte: on 
s'attend à ce que les événements rares créent une asymétrie dans la distribution 
de densité locale. Le problème est que le collapse sphérique ne tient pas compte 
des couplages non-locaux. Il ne saurait donc a priori donné une image fiable de 
ce qui se passe. L'objet des calculs qui suivent est d'explorer aussi précisément 
que possible ces couplages et leurs conséquences. 



6.1 Propriétés générales du développement perturbatif 

Je commence par donner des propriétés générales du développement perturbatif 
en reprenant des résultats obtenus par GorofF et al. (1984). Le point de départ 
est le système d'équation décrivant la dynamique du fluide en espace Eulérien 
dans l'approximation à un seul fluide. C'est d'ailleurs une hypothèse fondamen- 
tale qu'on ne sait pas lever. Je me place dans un espace Einstein-de Sitter pour 
faciliter la présentation (on peut le faire dans toute sa généralité mais c'est plus 
laborieux) . 

Je rappelle qu'on a, 

|p(x,t) + iv x .[p(x,t)u(x,t)] = (201) 
ot a 

d ' 1 1 

-u(x,t) + -u(x,f) + -(u(x,().V x )u(x,i) = — V x ^(x,t) (202) 

ot a a a 

VxV(x,t) -4xG\p(x,t) -p(<)]a 2 = 0. (203) 

A partir de là je vais prendre la divergence de la deuxième équation et utiliser 
l'équation de Poisson pour éliminer le potentiel. On a, 

Q 1 

-p(x,f) + -V x .[/>(x,i)u(x,t)] = (204) 
Ot a 

|- V x -u(x, t) + - V x -u(x, t)+ 

Ot a 

-^u iJ (x,t)u iJ (x,t) + -u(x,t).V x (V x .u(x,t)) = -~— 5(x,*X205) 

CL CL £t CL 

ij 

Enfin je vais introduire un nouveau gradient, V = Vx et le champ $ de telle 
manière 3 que u = V^. Enfin les dérivées par rapport au temps sont exprimées 
comme des dérivées par rapport à a. Notons qu'alors A<Ï>( X ) s'identifie avec 
9 qui est, on s'en souvient, une observable indépendante de la constante de 
Hubble. On obtient, 

d 

a— ôhc) + (1 + £(x))A$(x) + V£(x).V$(x) = (206) 
oa 



3 la partie rotationnelle du champ de vitesse correspondant à un mode décroissant, le champ 
de vitesse qu'on regarde est a priori potentiel. 
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d 1 
a— A$(x) + -A$(x) + V$(x) • VA$(x)+ 
oa 2 

Ç^^^W + l^x) - 0. (207) 

ij 

Si on écrit la densité (et la divergence) comme une série en fonction de la 
densité initiale, 

5(x) = 5>">(x), (208) 

n 

6>(x) = ^6>^(x), (209) 
et en introduisant les modes de Fourier de ce champ initial, 

*(*) = / j0r^^^ ( 21 °) 

alors on peut facilement montrer que, 

H 3 k 

j2^ S{kl) • ' ' j2^ S(kn) °" F " (kl ' ' ' ' ' k " } (2U) 

/rl 3 k H 3 k 
7^37^1) • • ■ J^pS(k n ) a" G n (k!, . . . , k„). (212) 

Les fonctions F n et G n sont des fonctions sans dimensions, homogènes, des 
vecteurs d'onde ki. Il existe des lois de recursivité liant ces fonctions entre elles. 
Par exemple, 

F l = 1 (213) 
d = -1 (214) 
5 lkj.k 2 lk!.k 2 2 (kx.k 2 ) 2 

"3 , lki.k 2 , lki.k 2 , 4(kx.k 2 ) 2 " 



G 2 = 



T 2 k-^ 2 7 fc-^ k<2 



(216) 



Dans le cas d'un Univers qui n'est pas EdS, cette propriété n'est plus valable 
au sens stricte. Mais si dans l'équation précédente on remplace les facteurs 
d'expansion a par D +1 les fonctions F n qui en résultent ne dépendent que très 
faiblement du temps (et donc des paramètres cosmologiques). Les fonctions G n 
sont essentiellement toutes proportionelles à /(fi). 

La conséquence de ces propriétés est que les développements perturbat- 
ifs qui se font formellement par rapport à la densité initiale, sont en fait des 
développements par rapport à la solution linéaire. 
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6.2 Un exemple: le calcul de F 2 et G 2 



A partir des équations couplées, on peut facilement écrire une relation récursive 
entre F\, G\, F 2 et G 2 . On a en effet, 



2 F 2 (k, k') + G 2 (k, k') + F l (k) Gi (k') 



1 + 



k - k' 



fc' 2 



-G 2 (k,k') + Gi(k) Gi(k') 



k • k' (k • k') 2 



fc' 2 



fc 2 fc' 2 



+ 2^(k,k') 



= 0(217) 
= 0(218) 



A partir de F\ — —G\ (qui s'obtient trivialement à partir de l'équation de 
continuité) on en déduit F 2 et G 2 . Ce calcul n'est valable que dans le cas 
Einstein-de Sitter. On voit facilement qu'on peut généraliser ce type d'approche 
à un ordre arbitraire. La dépendance avec les vecteur d'onde sera de plus en plus 
compliquée, mais fera toujours intervenir les 2 mêmes fonctions géométriques 
de base. 



6.3 Effet du couplage de mode: la skewness 

Une première conséquence de ce développement qu'il est facile de voir est 
l'apparition de couplages de mode. La densité à l'ordre 2 est une convoca- 
tion des modes linéaires. D'un point de vue observationnelle il est intéressant 
d'en déduire des propriétés statistiques a priori observables et un moyen de met- 
tre en évidence ces couplages de mode est d'examiner l'apparition de propriétés 
non-gaussiennes. Le premier moment non-trivial à apparaître est la skewness, 
moment d'ordre 3 de la densité. 

On cherche donc à calculer (5 3 ) et son terme dominant. On a 

(<5 3 ) = ((V 1) +<5 (2) + ...) 3 ). (219) 

Si on réordonne les termes qu'on obtient perturbativement on a, 

(<5 3 ) = ( (<5 (1) ) 3 ) + 3 ( (<5«) 2 ô^) + ... (220) 

Les termes suivants sont d'ordre plus grand en théorie des perturbations. 

Le premier terme de ce développement est identiquement nul pour des condi- 
tions initiales gaussiennes. Le terme suivant est donc a priori le terme dominant 
pour cette quantité. On a donc, 

(S 3 ) « 3((ô^y S^) (221) 

f d 3 kx f d 3 k 2 f d 3 k 3 f d 3 k 4 
" J (2tt) 3 / 2 ° J (2^) 3 / 2 a J (2tt) 3 / 2 a J (2^) 3 / 2 a X 

(J(ki) <5(k 2 ) <5(k 3 ) ô(k 4 )) F 2 (k 2 , ka) exp[i(k! + k 2 k 3 + k 4 ) • $!22) 
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Pour des conditions initiales gaussiennes il faut associer les modes de Fourier 
par paires. Si on associe k2 et k3 ensembles on a (à cause de F 2 ). Il reste 
deux termes qui donnent, 



<* 3 > 



6 



d 3 ki f d 3 k 



a 4 P(fci)P(fc 4 ) x 



(2tt) 3 J (2tt) 
5 l k!-k 4 1 ki ■ k 4 2 (ki ■ k 4 ) 2 
7 2 k\ ++ 2 kj + 7 k\k\ 



(223) 



En intégrant sur les angles entre ki et k 4 on obtient finalement (Pceblcs 1980), 

(S 3 ) = f<<5 2 ) 2 . (224) 
On voit qu'il apparaît un nombre pur, 34/7, et la tradition est de définir, 

_ (6*) 



Sa 



(225) 



Ce calcul cependant est académique parce que la densité a été prise ponctuelle- 
ment. En pratique les champs sont filtrés (que ce soit observationellement ou 
dans des expériences numériques) . Il faut alors tenir compte de cet effet dans le 
calcul de 53. La difficulté du calcul tient dans la complexité qu'il y a à intégrer 
la partie angulaire des vecteurs d'onde. Pour avoir la skewness de la densité 
locale filtrée, 6r, il faut en effet calculer, 

( 5 **> = 6|^|^| a 4 P(fci)P(fc 4 )^(fciiî)^(fc 4 iî)x 
1 ki • k 4 1 ki • k 4 2 (ki • k 4 ) 2 



2 k\ 



k>4 



7 



k- 2 h 2 



Wflki +k 4 |iî()226) 



La fonction W est la fonction de filtrage dans l'espace des phases. Elle dépend 
évidemment de la procédure de filtrage utilisée. Il se trouve que le résultat final 
peut prendre une forme simple pour un filtrage top-hat dans l'espace réel. Dans 
ce cas là, 



W(k) -- 

Pour cette fenêtre on a 
dfti2 



/ 



Ait 
jdïl 12 



47T 



W(|kl+k 2 | 

VK(|ki+k 2 



3tt J 3/2 (fc) 
2 fc 3 / 2 



(ki • k 2 ) 2 



1 + 



h 2 h 2 

ft l "-2 

ki • k 2 



— [sin(fc) — k cos(fc)] 
k 



W(k 1 )W(k 2 ) 



k'( 



3 

W(fci 



(227) 



(228) 



W(k 2 ) + -k 2 W'(k0(29) 



Ces propriétés s'obtiennent en utilisant le théorème de somme des fonctions de 
Bessel. On voit qu'on peut exprimer la fonction F 2 en utilisant les 2 expressions 
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10 

R (h- 1 Mpc) 
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Figure 16: Les paramètres S p pour 3 < p < 7. Comparaison des prédictions 
théoriques (Bernardeau 1994c) avec les résultats de simulations numériques 
(Baugh et al. 1994) 



polynomiales des relations précédentes. Finalement on obtient (Bernardeau 
1994b), 

S3 = f + i^>. (230, 
7 d log R 

Le résultat final va dépendre de la forme du spectre du puissance (essentielle- 
ment à l'échelle de filtrage). Pour un spectre en loi de puissance, 



P(K) ex k v 



(231) 



on a, 



34 

S 3 = y (n + 3) 



(232) 



L'étude des amas de galaxies permet d'avoir des indications sur la valeur de 
n. On trouve des valeurs de l'ordre de n « —1.5. Les comparaisons avec les 
simulations numériques ont montré que ce résultat perturbatif sur S3 était très 
robuste. Le domaine de validité de ce résultat est relativement grand comme 
on peut le voir sur la figure. 
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6.4 Dépendance avec les paramètres cosmologiques 



Une question très intéressante qui se pose alors est de savoir si ce paramètre S3 
qui quantifie l'émergence de propriétés non-gaussiennes dépend des paramètres 
cosmologiques. Les calculs analytiques ont montré que la dépendance attendue 
était très faibles. En particulier Bouchct et al. 1992 ont obtenu, 



S 3 = y + ^ -°- 03 -l)-(n + 3). 



(233) 



Ce résultat est basé sur le fait que la densité au 2ème ordre s'écrit maintenant, 



<5 (2) (x) = / A J A D% Sfr) 5(k 2 ) expftk! + k 2 ).x] x 



1 ^2 



(234) 



Ce résultat a été obtenu quand A 
A ^ 0. 



0. On trouve un résultat similaire quand 



6.5 Interprétation, effet du filtrage 

La skewness mesure la tendance du système à créer une distribution asymétrique 
avec des événements rares ayant un contraste de densité grand et positif (les 
proto-amas) et beaucoup d'événements de contraste de densité petit et négatif. 
Un moyen de quantifier cet effet est de considérer le développement de Edge- 
worth de la distribution de densité locale. Ce développement est valable pour 
des distributions proches d'une gaussienne (dans un sens qui sera préciser par 
la suite). Il s'écrit, 

^-tè? -*(-£) [ i+ f'*(i) + --- 

Ce développement utilise les polynômes de Hermite. Ici 

H 3 (y) = v 3 -Zv. (236) 

Sa forme générale peut être vue comme un développement formel (et par in- 
dentification on peut relier les coefficients des termes du développement avec 
les moments, Juszkiewicz et al. 1995) ou à partir de la forme formelle de la 
distribution de densité locale obtenue comme transformée de Laplacc inverse de 
la fonction génératrice des cumulants (Bernardeau & Kofman 1995). 

Pour comprendre la dépendance de ce paramètre avec la forme du spectre il 
est très intéressant d'examiner en détail la nature du terme qui est rajouté avec 
les effets de filtrage. 

Un bon moyen d'appréhender ce problème est de regarder ce qui se passe 
dans l'espace Lagrangien. En effet si on calcule on obtient, 

J(2) = S (07^ / ^n aH{ ^ ^ -P[Kk 1 + k 2 ).q] x 



(235) 
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Figure 17: La skewness est une mesure de l'asymétrie de la distribution de den- 
sité locale. Elle émerge parce que les régions sous-denses évoluent moins vite 
que les régions sur-denses dès que les non-linéarités commencent à jouer un rôle. 
La dépendance de la skewness avec la forme du spectre vient de la correspon- 
dance entre l'espace Lagrangien qui donne la taille initiale de la perturbation 
et l'espace Eulérien qui est sensible à la taille finale. A une échelle de filtrage 
R donnée, les régions sur-denses proviennent de l'effondrement de régions qui 
étaient initialement plus grandes, alors que les régions sous-denses détectées vi- 
ennent de régions initialement plus petites. Du coup si on part d'un spectre 
de puissance avec beaucoup plus de fluctuations à petite échelle qu'à grande 
échelle, on s'attend a ce que la skewness soit plus petite. 
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2 
7 



(k!.k 2 ) 2 

h 2 h 2 



(237) 



ce qui donne pour la densité après que le Jacobien (c'est à dire le volume) ait 
été filtré à une échelle R donné, 



f(2) 



(ki.k 2 ) 2 



VK(fci i?) VK(fc 2 J2) - -W(|ki + k 2 1 iî) 



h 2 k 2 



Il en résulte que le paramètre 53 devient, 

^Lag. _ 



34 



(238) 



(239) 



Le fait que l'on ne trouve pas le même résultat ne doit pas surprendre. Dans 
ce cas on vient de faire un filtrage à une échelle de masse donnée. La différence 
entre les 2 c'est qu'un objet de masse donnée a a priori une taille d'autant 
plus petite que cette masse était grande. Un filtrage à une échelle Eulérienne 
donnée mélange donc différentes échelles de masse initiales. Or pour un spectre 
hiérarchique les fluctuations de masse sont d'autant plus faibles que la masse est 
grande. On aura donc une asymétrie moins importante que celle qu'on aurait 
pu escompter. 



7 La hiérarchie des corrélations en régime quasi- 
linéaire 

Peut-on faire mieux que le calcul de S3? Le calcul des fonctions homogènes 
F n ou G n est d'autant plus difficile que l'ordre est élevé, et les intégrales qu'il 
faut calculer deviennent elles aussi très fastidieuses. En pratique on ne peut 
guère aller au delà du moment d'ordre 4 par un calcul direct. Fry (1984) avait 
explorer ce cas mais sans prendre en compte les effets de filtrage. Bernardeau 
(1994b), Lokas et al. (1995) ont fait ces calculs directs pour S 4 dans le cas 
respectivement d'un filtre top-hat et d'un filtre Gaussien. 
Formellement le cumulant d'ordre 4 est donné par, 

(S 4 ) c = (<5 4 >-3<5 2 ) 2 (240) 

= 12<(*W) a (^) 2 ) c + 4((^)) 3 ^)) c . 

Dans ces équations le fait de prendre la partie connexe devient essentiel. Les 
termes supplémentaires qui apparaissent sont des corrections à la variance. 
Ils doivent donc être enlevés (ils s'annulent naturellement quand on fait la 
différence). 
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Une conséquence est que, 



(241) 



et on peut définir S 4 par 



2\3 



^ - <<5 4 ) c /<* 2 ) 



(242) 



L'équation précédente permet de calculer la partie dominante de £4 en régime 
quasi-linéaire. D'une manière générale on définit 



S n = (ô n )J{5 2 ) 



2\(«-l) 



(243) 



Toutes ces quantités sont finies à grande échelle pour des conditions initiales 
Gaussiennes. L'objet de cette section est d'en faire le calcul. 



7.1 La fonction génératrice des cumulants de la densité 
locale 



On peut déjà remarquer que 



5 4 = 12^| + 4^3, 



(244) 



avec 



v 2 



(6™ \SM] )/< \SM] ) 



tf 1 



>e/( 



5^ 



) 3 - 



(245) 
(246) 



d'une manière générale tous les S n peuvent s'exprimer en fonction des seuls 
quantités, v p , avec 



(247) 



Pour le voir on peut faire une représentation graphique des termes qui con- 
tribuent à S n . Dans chaque 6^ il y a un produit de p variables aléatoires 
gaussiennes <5(k). Je vais représenter chacune de ces variables par des points. 
Quand on doit calculer une valeur moyenne, l'application du théorème de Wick 
fait que tous ces points doivent être associés par paires. Les quantités 
ressemblent donc à ce qui est représenté sur la figure |î^. 

Pour calculer S n il faut calculer le terme dominant de la partie connexe de 
produits de la forme ô^ 1 ^ . . .ô( Pn >. Pour que le terme calculé contribue à la 
partie connexe il faut que tous les facteurs soient liés par au moins une patte. 
Si ce n'est pas le cas on obtient un terme qui entre dans les termes correctifs au 
cumulants d'ordre inférieur. Par ailleurs comme on cherche le terme dominant 



56 




etc.. 



Figure 18: Représentation diagramatique des S p . Chaque point représente un 
facteur <5(k). 




Figure 19: Expression des graphes les plus simples. Chaque ligne porte un 
facteur a 2 , Les vertex obtenus par moyenne angulaire sur les vecteur d'onde 
portent v p . 




Figure 20: Un exemple de graphe contribuant à S$. 
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Figure 21: Représentaion graphique de la relation fl252| ), r est la fonction 
génératrice des diagrammes à une patte externe. 

le nombre de lignes utilisées doit être minimal et le nombre de lignes nécessaire 
pour relier n points étant n-1 on aura, 

S n = J2 (ô M ...S^) c /{[ô^] 2 ) n -\ (248) 

graphes, Pi— 2(n— 1) 

Un exemple d'un tel graphes pour S§ est donné sur la figure 

Il est utile de remarquer que tous ces diagrammes sont des arbres, qu'ils 
ne contiennent pas de boucles. On peut alors intégrer sur les vecteurs d'onde 
de proche en proche. Chaque ligne fait apparaître un facteur a 2 et chaque 
vertex fait apparaître un facteur v v tel que je les ai définnis précédemment. 



Graphiquement cela donne ce qui est représenté sur la figure 19. 

Ce qu'on aimerait pouvoir calculer c'est l'ensemble des coefficients S n) et 
donc plus spécifiquement leur fonction génératrice f{y), 

V (y) = jr / -S n i -^-. (249) 

Si on définit la fonction génératrice des vertex Gs(t) par, 

g s (T)=J2"n ( -^, (250) 
P =i 

il est possible de montrer ip et Qs sont liés par le système, 

v(v) = y + yG5[T(y)] + \r 2 (y), (251) 

(y) = -vG'Mv)]- ( 252 ) 



T 



Pour démontrer ce système de relations il faut faire intervenir les transformées 
de Legendre des fonctions génératrices. Ces techniques ont été développées 
initialement dans le domaine des polymères (voir Des Cloiseaux & Jannink 
1987). Dans un contexte de cosmologie elles ont été employées par Schaeffer 
(1984) et présentées en détails par Bernardeau & Schaeffer (1992). S'il n'est 
pas possible de démontrer cette relation en quelques lignes on peut au moins se 
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convaincre de la validité de l'équation ( |252| ). En effet r est la fonction génératrice 
des graphes à une patte externe. Sur la figure ( pï[ ) on voit que r doit être solution 
d'une équation implicite qui fait intervenir la fonction génératrice des vertex. 

Finalement la difficulté technique est dans le calcul de la fonction génératrice 
des vertex Gs(t). Il faut avoir bien à l'esprit que les vertex v p sont simplement 
les moyennes géométriques des fonctions F p . Pour faire les calculs je vais aussi 
avoir besoin des moyennes géométriques des fonction G p que je note /i p , soit, 



■angles 



G p (ki, . . . ,k p ) 



■angles 



(253) 
(254) 



et je défini Çg{r) la fonction génératrice des vertex \i p . Pour obtenir Qs et Ge il 
faut partir des équations (206-207) et les ré-écrire dans l'espace de Fourier, 



n-l 



nF„(ki 



■ ! k n ) 



•Gn(kl 



• i k n ) 



• kp) (x n . p(kp k n ) x 



0=1 



i + * — ^ — 



\£~ii= p +l 



= (255) 



Gn(ki 



.,k rs ) + -iî , n (k 1 ,. 



■ ,K) 



n-l 



^ ] Gp (ki , . . . , kp) G n _p(kp_|_i , . . . , k n ) x 
P =i 

(Ei=i ki) ■ fSi=p+i k^ [(Si=i ki) • (j2i- p+ i ki 



Si=p+i k 



Cr=iki) (Z)i=p+ik 



= 0. (256) 



Quand on prend les moyennes angulaires les produits scalaires simples tombent 
tandis que les produits carrés donnent 1/3. On obtient donc, 



n v n + n n + y v p fi n - 



p=l 



pir, 



1 " 



0: 



= 0. 



p=\ 



En termes de fonctions génératrices cela implique que, 

T^-g s (T) + g e (r) [i + Gs(t)} 

dr 



0: 
0. 



(257) 
(258) 

(259) 
(260) 
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On obtient un système d'équations différentielles couplées. Après quelques ma- 
nipulations on peut montrer que l'équation d'évolution du Qs(t) se réduit à celle 
de l'effondrement spherique où G s (t) est le contraste de densité de la structure 
en effondrement et r est la sur-densité linéaire. Dès lors il facile de calculer les 
v n à un ordre arbitraire, 

34 , s 

^2 = (261) 

^3 = — ; (262) 
A 189' y ' 

446440 ,„„„, 

^ = ■; ( 263 ) 



43659 



Les valeurs des S n en découlent tout de suite, 



(264) 



34 

S 3 = -■ (265) 

Si = l^" 45 - 89; (266) 

S 5 = g^Sf ~ 656.3; (267) 
305 613 v ' 

(268) 

C'est intéressant mais il reste que ces résultats ne tiennent pas compte de l'effet 
de filtrage. 

7.2 La fonction génératrice des cumulants de la densité 
filtrée 

7.2.1 Propriétés géométriques dans l'espace Lagrangien 

La quantité de base en espace Lagrangien est le Jacobien. Je peux toujours 
imaginer un développement perturbatif de J, 

j( q ) = J (1) (q) + J (2) (q) + --- (269) 

A un ordre donné on aura, 

/d 3 ki d 3 k 
(2^ ' ' • (2^ ° P Jp(kl ' ' ' ' ' kp) cx P^-( kl + • ' • + k ^ ( 27 °) 

La propriété remarquable que l'on peut montrer est que (Bernardeau 1994c) 

j P = J P (k 1 ,...,k p )T^(|k 1 + ... + k p |iî) anslos (271) 

-angles 



J p (ki,...,kp) W(feii2)...W(fepiî). (272) 
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La démonstration complète est assez ardue. Elle repose sur les remarques 
suivantes. Le Jacobien s'exprime essentiellement à partir de deux fonctions 
géométriques dans une construction récursive. En effet 



9x 



= 1 



Vq.4< 



■*) a -E*w*A< 



(V q .ï>f - 3V q J2 (| '- +2£*y*jAi 



tjk 



(273) 



On a déjà vu que le filtre top-hat en quelque sorte commute avec 1— (ki.ka) 2 /fci/A^. 
On peut aussi montrer la propriété suivante, 



/ dftidft 2 dft 3 W (|ki + k 2 + ks| fl) 



1 - 



ki.k2 



_ ( k 2 .k 3 \ _ / k 3 .ki \ i o ki.k2 k 2 . 

(4tt) 3 § W{k x R) W(k 2 R) W(k 3 R). 



k3 k3.ki 




(274) 



Là encore, tout se passe comme si on prenait le résultat obtenu pour le jacobien 
ponctuel calculé à partir du champ de densité initial tronqué. 

Il est relativement facile de voir que la fonction J p n'est construite récursivement 
qu'avec de telles fonctions géométriques. La propriété ( |272| ) s'obtient alors par 
applications successives de ces propriétés géométriques. 

Ce qu'on vient d'obtenir généralise ce que l'on avait vu pour S3 où en effet 
5*3 dans l'espace Lagrangien n'était pas affecté par le filtrage (pour un filtre 
top-hat en tout cas). 



7.2.2 Le passage Lagrangien-Eulérien 

L'idée est maintenant de trouver un truc pour récupérer le terme dominant des 
S p à moindre frais à partir des S p de l'espace Lagrangien. 

La hiérarchie précédente me donne des informations sur la fonction génératrice 
des cumulants de la distribution de volume à échelle de masse fixe. On peut 
alors faire la remarque suivante: la probabilité qu'une masse M occupe un vol- 
ume supérieur à V est aussi la probabilité qu'un volume V contienne une masse 
inférieure à M, Il suffit pour cela de considérer des sphères concentriques autour 
d'un point x donné^. Or une de ces probabilités s'exprime en espace Lagrang- 
ien, l'autre en espace Eulerien. A travers cette proposition on peut alors trouver 
le comportement dominant de la fonction génératrice des cumulants <p(y) dans 

4 La proposition n'est cependant rigoureuse que pour des probabilités dites centrées 
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la limite des petites variances à partir de celle obtenue pour l'espace lagrangicn 
qui s'identifie à celle obtenu sans tenir compte des effets de filtrage. 
Rappelons d'abord que 

+ioc r , ys- dS _ (275) 



P(ô)dô = J — ,ex P 



Il est alors intéressant de faire un calcul succinct par la méthode du col qui est 
a priori valable dans la limite des petites variances. On voit que la position du 
point col est donnée par l'équation 

^(V)=S, (276) 

et par ailleurs on a 

^<P(v)=G S (r), (277) 

quand r est toujours donnée par la même relation implicite avec y. Autour de 
ce point col la partie sous l'exponentielle est simplement donnée par 



exp(-^). (278) 



La position du point col est donc essentiellement donnée par un simple change- 
mement de variable passant de la densité linéaire r au contraste densité non- 
linéaire S. 

La différence entre une description Lagrangienne est que a est pris dans 
un cas à une échelle de masse donnée, dans l'autre à une échelle physique R 
donnée. Comme je cherche la partie dominante de ip dans la limite a — > il est 
parfaitement légitime de s'appuyer sur cette approximation. 

Pour rendre le calcul plus simple, je vais préciser un peu mes notations. Je 
note ip L la fonction génératrice des cumulants dans l'espace Lagrangien et Gs 
la fonction ghératrice des vertex. Dans l'espace Eulérien j'utilise l'exposant E. 
Il est toujours possible de supposer qu'il existe une fonction Gf associée à <p E 
(même si il n'y a pas de représentation diagrammatique à partir de vertex qui 
lui soit associée). La variable r est de même prise avec l'indice L ou E selon 
qu'elle concerne l'espace Lagragien ou Eulcrien. 

Par identification des termes sous l'exponentielle (les seuls qui comptent dans 
la limite considérée), on a 

(279) 



2<j 2 (R) 2a 2 [(1 + S) 1 / 3 R] ' 

Le contraste de densité est un paramètre donné a priori. Les variables te et tl 
dépendent de 6 à travers les équations du col, 

QHtl) = ô, (280) 
Gf(r E ) = ô. (281) 
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Dans le cas Lagrangien a (fonction de l'échelle à priori connue) est pris à l'échelle 
de masse correspondant à la surdensité ô considérée. La fonction Gs( t l) est 
connue. C'est elle qui décrit la dynamique de l'effondrement sphérique. 

A partir ce ces équations je peux éliminer tl pour obtenir une équation 
implicite entre Qf et te, 

5s (te) = Qs y ^ r E j . (282) 

La fonction génératrice des cumulants (p E (y) est alors construite à partir de 
G g (te) de la même manière que tp L (y) l'était à partir de Gs( t l)- 

En développant cette fonction autour de y = on a alors les expressions 
explicites des prmières valeurs des S p . Ils s'expriment en fonction des dérivées 
logarithmiques successives de la variance, 

_ d'iog a 2 (R) 



dlog p i? 



et alors, 



5 3 - y +7i, (284) 

60712 62 7l 7 7l 2 2 72 

54 ~ T323" + ^~ + ^~ + ~' (285) 



S: 



200575880 1847200 7l 6940 V 235 7l 



305613 3969 63 27 

1490 72 50 7i 72 10 73 



(286) 



63 9 27 ' 

= 12650 + 12330 71 + 4512 7l 2 + 734.0 7i 3 + 44.81 71 4 
+775.8 72 + 375.9 71 72 + 45.56 71 2 72 + 3.889 7 2 2 
+20.05 73 + 4.815 7i 73 + O.I85274, (287) 
S 7 = 307810 + 383000 71 + 190700 7 i 2 + 47460 71 3 

+5914 7l 4 + 294.8 71 5 + 27340 72 + 20300 71 72 

+5026 71 2 72 + 414.8 71 3 72 + 358.1 72 2 + 88.15 71 7a 2 

+902.6 73 + 443.3 71 73 + 54.44 71 2 73 + 7.778 72 73 

+ 14.20 74 + 3.4577i 74 + 0.08642 7s, (288) 



La figure ( |Ï6[ ) montre que ces coefficients sont en parfait accord avec des 
résultats de simulations numériques (ici pour un spectre CDM) dès lors que 
l'échelle de filtrage est plus grande que 10/i -1 Mpc. Il est assez remarquable 
de voir que le domaine de validité de ces résultats ne se détériore pas quand on 
augmente l'ordre des cumulants. 
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Table 1: Valeur des paramètres de la singularité (|293| ) pour différentes valeurs 



de l'indice spectral n. 






n 


Sa 


Vs 


a s 




-3 


0.656 


-0.184 


1.84 


-0.030 


-2.5 


0.804 


-0.213 


2.21 


-0.041 


-2 


1.034 


-0.253 


2.81 


-0.058 


-1.5 


1.44 


-0.310 


3.93 


-0.093 


-1 


2.344 


-0.401 


6.68 


-0.172 


-0.5 


5.632 


-0.574 


18.94 


-0.434 



7.3 Application, construction de la distribution de prob- 
abilité de densité locale 

Je vais reprendre la relation entre la fonction de distribution de probabilité et la 
fonction génératrice ip (y). Pour pouvoir utiliser cette relation il faut faire une 
hypothèse supplémentaire non triviale. En effet ip(y) est a priori une fonction 
de ex à travers la dépendance des S p avec a. Je vais alors supposer que la limite, 



Vil), c) ->• f(y) quand er -> 0, 



(289) 



est uniforme. Il n'y a pas de preuve de cette propriété, cependant elle est 
fortement suggérée par les résultats numériques sur les S p . 

Quand la variance est suffisemment faible il est alors légitime de supposer 

que 

<p(y) , yà~ 



f +i °° ày 

pms = L i0O w exp 



dô. 



(290) 



Je vais maintenant considérer sérieusement cette équation et rechercher les 
différentes formes de p(ô) qui en résultent (voir Balian & Schaeffer 1989 pour 
une description plus exhaustive des résultats). Profitant de l'approximation 
a <C on peut appliquer la méthode du col, et on obtient 



p(ô)dô 
5s(r) 



dô 



l-^'(r)M(r) 



-,1/2 



27TCT 2 



exp 



2a 2 ) 



= S. 



(291) 
(292) 



Cette solution est valable quand S < ô c où S c est la valeur du contraste de 
densité qui annule 1 — tÇ' s '(t)/G' s {t). Quand 5 est plus grand que ô c on ne 
peut plus faire l'approximation du point col. La forme de p(5) est alors plutôt 
déterminée par le comportement de ip(y) près de sa singularité, 



tp(y) - tp s ~ -a s (y - y s ) 



3/2 



et on a 



P (ô)dô 



3a s a 



(1 + Ô)- 5 / 2 exp [-|y s |<ï/<7 2 + \^ s \/a 2 } dô. 



(293) 



(294) 
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Figure 22: Comparaisons des résulats de simulations numériques N corps avec 
les prédictions de la théorie perturbative. 



On voit que la forme du cut-off est très différente d'une forme gaussienne. Il faut 



bien avoir à l'esprit que ce résultat repose en grande partie sur l'hypothèse (289), 
et en particulier que la position de la singularité doit rester à une distance finie 
de l'origine quand a est finie. La table (|ï|) donne les paramètres de la singularité 
pour différentes valeurs de l'indice spectral. Je rappelle que n sa —1.5 est une 
valeur suggérée par les comptages d'amas. 

Numériquement il est toujours possible d'intégrer sans faire explicitement 
l'approximation du col. Les formes obtenues sont en bon accord avec les résultats 
de simulations numériques ./V-corps. Sur la figure ( p2|) on a des fonctions de 
distribution pour différentes échelles de filtrage, donnant différentes valeurs de 
la variance. La forme prédite pour la distribution (calculée à partir de la donnée 
de la variance et de l'indice spectral local) est en remarquable accord avec les 
résultats numériques. 



7.4 Les termes sous-dominants du développement pertur- 
batif 

Les quantités que l'on vient de calculer sont les termes dominants (contributions 
en arbre) d'un développement plus général. On peut toujours se demander si les 
corrections en boucles permettent d'étendre le domaine de validité de la théorie 
des perturbations. Par exemple pour les termes correctifs à la variance on a 
deux diagrammes supplémentaires (voir figure ^3j), pour la skewness on en a 
quatre (figure [?4|) . 



G5 




Figure 23: Diagrammes à une boucle contribuant à la variance 




Figure 24: Diagrammes à une boucle contribuant à S3. 



C'est un domaine qui a été relativement peu exploré. On sait donc peu de 
chose sur ces termes correctifs. La raison principale en est sans doute que les 
calculs sont horriblement compliqués. On ne peut plus simplement moyenner 
sur les angles, et les effets de filtrage rendent les calculs encore plus inextricables. 

Les résultats obtenus sur la variance et la skewness montrent que (Scocci- 
marro & Frieman 1996, Scoccimarro 1997, Scoccimarro et al. 1998) 

• Les termes correctifs ne convergent que si l'indice spectral n est compris 
entre -3 et -1. A grande échelle ça ne converge donc pas. Ce sont les 
'petites boucles' qui induisent ces divergences. 

• Quand on a convergence, on a une bonne description de la transition vers 
le nonlinéaire (voir figure p5|) . 

En pratique on a un cut-off naturel parce que n — > —3 à petite échelle. Il n'en 
reste pas moins qu'on ne sait pas si ce cut-off doit être explitement introduit 
pour régulariser les résultats ou si une resommation (au moins partielle) des 
boucles donnerait une contribution finie (mais éventuellement avec une autre 
dépendance en a). 
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Figure 25: Résultats de la théorie des perturbations à une boucle pour la vari- 
ance et le bi-spectrum pour le cas n = — 2. On voit que les termes en boucles 
reproduisent correctement les résultats numériques, (figure extraite de Scocci- 
marro et al. 98) 
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Figure 26: La skewness mesurée dans le catalogue IRAS (carrés) comparée à des 
mesures faites pour les galaxies optiques (triangles). La ligne continue à gauche 
est la prédiction de la théorie des perturbations (figure extraite de Bouchet et 
al. 1993). 

8 Application aux propriétés statistiques de quan- 
tités observables 

Dans cette partie je vais illustrer les résultats des théories des perturbations 
sur des quantités observables, catalogues de galaxies tri-dimensionnel ou bi- 
dimcnsionncl, champ de distorsion gravitationnelle. 

8.1 Les catalogues de galaxies 

A priori on peut toujours supposer que les fluctuations du nombre local de 
galaxies sont une mesure des fluctuations de densité cosmique. En général on 
suppose que les galaxies forment une représentation discrète résultant d'un pro- 
cessus Poissonnien d'un champ de densité continu sous-jascent. (Ce qui revient 
à dire que les mesures sont affectées d'un bruit Poissonnien). 

Dans les catalogues tri-dimensionnels il faut tenir compte d'un effet ignoré 
jusqu'à présent qui est dû au fait que les positions des galaxies sont estimées à 
partir de la vitesse de recession qui combine la distance réelle (mouvement de 
Hubble) et la vitesse propre le long de la ligne de visée. C'est ce qu'on appelle la 
position dans l'espace des redshifts. Il faut en principe en tenir compte dans le 
calcul des S p . Hivon et al. (1995) ont montré que cette correction ne changeait 
que très peu la valeur de S3 à grande échelle (enfin tant qu'on peut négliger les 
'doigts de Dieu', voir figure |27| ). 

Des mesures de ces paramètres ont été faites dans tous les grands catalogues 
de galaxies à notre disposition par Bouchet et al. (1992, 1993 pour le survey 
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log 10 ° Z 

Figure 27: La skewness calculée dans l'espace réelle (ligne pointillée) et calculée 
dans l'espace des redshifts (ligne continue). Ces prédictions sont comparées à 
des résultats de simulations numériques, espace réel (symboles pleins) espace 
des redshift (symboles vides). La cas présenté est n = — 1 pour un filtre top-hat 
(figure extraite de Hivon et al. 1995). 

de galaxies infrarouges IRAS), Gaztanaga (1992, 1993, et 1994 pour les surveys 
optiques CfA, SSRS et APM). 

Toutes ces mesures montrent que les comportements hiérarchiques atten- 
dus sont bien vérifiés. C'est en soit un résultat extrêmement important. Cela 
implique que les conditions initiales pour la formation des grandes structures 
étaient sinon Gaussiennes du moins proche d'être Gaussiennes. Par exemple 
les données excluent a priori les modèles de textures (Gaztanaga & Mâhôncn 
1996, voir figure |2^). En effet on s'attend à ce que la hiérarchie soit préservée 
même si les galaxies sont biaisées. (du moins tant que les intéractions non- 
gravitationnelles ayant donné naissance aux galaxies sont locales). 

Toutefois les valeurs des S p mesurés ne sont pas en accord avec les prédictions 
théoriques. C'est à ce niveau là que se pose le problème des biais, c'est à ce jour 
l'interprétation la plus raisonnable pour expliquer les différences persistantes qui 
subsistent. 

8.2 Les catalogues bi-dimensionnels 

Les données les plus riches actuellement sont dans les catalogues angulaires. Ils 
donnent la position angulaire sur la voûte céleste de plus d'un million d'objets. 
De ce fait il est possible de mesurer les paramètres s p jusqu'à de très grands 
ordres (voir figure ||^) . 

On peut tout aussi bien se poser la question du calcul des probabilités de 
comptages dans ce type de catalogues. Cela correspond à un filtrage qui n'est 
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Figure 28: La skewness pour des modèles à conditions initiales Gaussi- 
ennes (symboles vides) compare à des modèles avec conditions initiales non- 
Gaussiennes de type textures (symboles pleins) pour différentes valeurs de la 
variance. Figure extraite de Gaztanaga & Màhônen (1996). 
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Figure 29: Les paramètres s p mesurés dans le catalogue angulaire APM 
(Gaztanaga 1993). 
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plus à symétrie sphérique: la fenêtre de filtrage a maintenant la forme d'un cône 
dans l'espace réel. Plus précisément la densité projetée w(j) dans la direction 
7 vaut, 

/>oo 

10(7)=/ r 2 drF{r)p(l + 6(r,j)), (295) 
Jo 

et la densité filtrée à une échelle angulaire 9 est 

Ml)= [ d V / r 2 drF(r)p(l + «5(r,7 + 7')). (296) 
J\Y\<e Ja 

La fonction F est directement liée au processus de sélection des galaxies du 
catalogues (limite en magnitude, etc.). Typiquement pour des sélections en 
magnitude apparente en luminosité optique on a, 

F(r) ~ cxp(-r 2 /^ 2 ), (297) 

où D est la profondeur du catalogue. Notons que jusqu'à présent les catalogues 
sont peu profonds d'un point de vue cosmologique et donc on peut supposer que 
la distance lumineuse est une simple distance Euclidienne. On peut toujours 
normaliser F de telle manière que 

/ r 2 drF(r)p = 1. (298) 
Jo 

Je cherche maintenant à exprimer les différents moments de la fonction de 
probabilité de densité angulaire locale. La fonction de corrélation à deux points 
est donnée par, 

M71M72)) = J rîdnF(n) j r 2 dr 2 F(r 2 ) J -0 P(k) x 

exp [ik r (n - r 2 ) + ikj_.(ri7i - r 2 7 2 )] , (299) 

où k est décomposé en sa partie radiale k r et sa partie ortho-radiale kj^ . L'approximation 
des petits angles consiste à remarquer que la distance angulaire I71— 72I est faible 
alors k r va être très petit devant la norme de kj_, du coup, 

P(k) « P(jfex). (300) 

Il en résulte que 

(^(71)^(72)) = j r 4 drF 2 {r) j Ajl p(fc) exp [ik ± .( 7 i - 72 )r] . (301) 

En effet l'intégrale sur k r introduit une fonction delta en r\— r 2 . L'approximation 
des petits angles signifie simplement que les paires qui contribuent le plus à la 
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variance de la densité projetée viennent de points qui étaient à la même dis- 
tance de l'observateur. Ce résultat correspond à l'approximation dite de Lim- 
ber (1954) sur la focntion de corrélation angulaire. Cette approximation n'est 
d'ailleurs valable que si la fonction de corrélation est suffisamment pentue. 

Pour avoir la variance de la densité filtrée, il suffit de prendre la TF de la 
fenêtre donc, dans l'espace de Fourrier, multiplier par la fenêtre^] W(k) et on 
obtient, 

(«$) = J r 4 àrF 2 {r) J AJl p(fc) W 2 (k 9 r). (302) 

Le calculs de la skewness peut se faire de manière similaire et on peut utiliser 
pour un filtrage angulaire top-hat de propriétés similaires à celles utilisées pour 
le cas 3D. 

Il en résulte que, 



S 3 



36 

y 



fr 6 drF 3 (r)f^P(k) W 2 (k6r) 



d 2 k 

(2tt) 2 



P(k')W 2 (k'6r) 



2 logé» J (2tt) 



P(k') W 2 {k'9r) 



(.303) 



Cette expression se simplifie quelque peu quand on suppose qu'on a une loi 
de puissance pour le spectre P(k) (P(k) ~ k n ), 



S3 



J r 6 àrF 3 (r) 



„-2(n+2) 



[JVdrF 2 (r) r-(™+ 2 )] 



36 

y 



-(n + 2) 



(304) 



La hiérarchie complète des s p se calcule de manière analogue au cas 3D. On 
a simplement en plus un facteur géométrique lié aux effets de projection. Les 
coefficients purement dynamiques sont en fait identiques à ceux qu'on aurait 
obtenus si on était parti d'une dynamique bidimensionnelle. 

Encore une fois, les comparaisons avec les catalogues de galaxies montrent 
un désaccord avec les prédictions de la théorie des perturbations, probablement 
dû aux effets de biais. 



8.3 Les effets de lentille gravitationnelle 

Une motivation pour aller chercher de nouveaux traceurs du champ de densité 
est justement de s'affranchir des effets de biais. Un moyen de tracer les champs 
de densité locaux est d'utiliser les effets de distorsion gravitationnelle. Pour une 
description général de ces effets voir Sachs (1961) et Schneider, Ehlers & Falco 
(1992). 

5 A deux dimensions la fenêtre top-hat de l'espace réel s'écrit en espace de Fourier, W(k) = 
2Ji(fc)/fc 
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Figure 30: Les paramètres S3 et S4 obtenus par la théorie des perturba- 
tions (lignes continues) comparés à des résultats de simulations numériques 
(Gaztanaga & Bernardeau 1998). 
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8.3.1 Les lentilles gravitationnelles 

Je vais commencer par rappeler des concepts simples des effets de lentilles. Une 
lentille induit une déflection da des rayons lumineux pau unité de distance ds 
égale à, 

da/ds = -2V ± </> (305) 

où <j) est le potentiel de la lentille et le gradient est pris dans le plan orthogonal à 
la trajectoire des photons. Quand on regarde des objets lointains ces déflections 
induisent des déplacements des positions apparentes des objets d'arrière plan. 
Le déplacement apparent dépend des différentes longueurs qui interviennent 
dans le banc optique. 11 en résulte que la position apparente d'un objet j 1 
(dans le plan dit image) est liée à sa position réelle dans le plan source j s par, 

7 J =7 S + £(7 J ) (306) 



avec 



^ I) = - 2 Do S L Dol W ^- i;J,,,7) 



où Dls, Dos et Dç>l sonr respectivement les distances angulaires entre la 
lentille et les sources, l'observateur et les sources et l'observateur et la lentille. 
D'autre part V 7 ^ est la dérivée angulaire du potentiel projeté $. 

Cependant le champ de déplacement n'est pas une observable directe (sauf 
dans les cas d'images multiples). Un objet (faiblement) étendu comme une 
galaxie va être non seulement déplacé mais aussi déformé. La matrice de 
déformation est donnée par la matrice 2x2 des dérivées secondes du champ 
de déplacement, 

Dt q 

A = Ç iti = S t3 - ^— (308) 
On voit que la matrice A est une matrice symétrique. On l'écrit en général, 

A=( l - K -^ ). (309) 

La trace de la matrice donne la convergence k, les autres composantes donnent 
le cisaillement. On va avoir des lignes critiques (et des images multiples) quand 

K — > 1. 

8.3.2 Les observables 

Pour chaque galaxie on peut calculer une matrice de forme, 5, à partir de la 
distribution de lumière de cette galaxie. Plus précisément on définit, 



(310) 

/ à 2 ee 2 x{9)/2 
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En se souvenant que la densité de luminosité d'un objet n'est pas affectée par 
les effets de lentille on peut relier la matrice de forme d'une galaxie d'arrière- 
plan telle qu'elle peut être mesurée dans le plan Image S 1 à celle qu'on aurait 
obtenue en l'absence de lentilles, 

det(^) v ; 

Chaque galaxie d'arrière plan permet donc une estimation de la matrice de 
distorsion, mais normalisée de telle sorte que son déterminant vaut 1. Du coup 
les matrices de forme donnent accès aux quantités suivantes, 

Cependant une galaxie d'arrière plan seule ne permet pas de faire une mesure 
(sauf cas d'arcs critiques ou quasi-critiques). Il faut moyenner un certain nom- 
bre de galaxies d'arrière-plan en supposant que les ellipticités intrinsèques des 
galaxies sont distribuées de manière aléatoire et indépendante. 



8.3.3 La précision des mesures possibles 

Le bruit intrinsèque des mesures est un bruit de type Poissonnien qui dépend 
de l'ellipticité intrinsèque des galaxies d'arrière-plan, donc 



Bruit ~ 




Typiquement on détecte environ 50 gal/arcmin 2 . À une échelle de l'ordre de 10 
aremin on a donc un bruit de mesure intrinsèque de l'ordre de 5 10~ 3 . Autrement 
dit on est en mesure de détecter des signaux cosmologiques à cette échelle dès 
qu'ils atteignent l'ordre du %. 

Un amas de galaxies génère des effets de l'ordre de quelques % à sa périphérie. 
Les effets de distorsion qu'il induit, même en dehors des régions critique, sont 
donc parfaitement détectables. En fait même les grandes structures de l'Univers 
en générale sont parfaitement détectables comme l'ont démontré Blandford et 
al. (1991), Miralda-Escudé (1991), Kaiser (1992) et comme on va le voir dans 
la suite. 



8.3.4 la relation entre cisaillement et convergence 

En tirant profit du fait que les éléments de la matrice A sont des dérivées 
secondes d'un même champ scalaire <j) on peut trouver une relation non-locale 
entre la convergence locale et le champ de distorsion. Cette relation explicitée 
par Kaiser (1995) s'écrit, 
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VIog(l-«)= - - À * • r (314) 
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Figure 31: Exemple de reconstruction de carte de convergence. La carte de 
gauche est une carte simulée de 5x5 degrés. La carte de droite est la carte 
reconstruite en introduisant du bruit de mesure pour chaque pixel. On voit à 
petite échelle le bruit Poissonnien des mesures mais à grande échelle on retrouve 
bien les grandes structures (figure extraite de Van Waerbeke et al. 1998). 

De cette relation il résulte que le champ de convergence locale est a priori 
observable. En fait, en pratique on ne fait pas une intégration directe de cette 
relation mais on recherche le potentiel <p qui reproduit au mieux le champ de 
distorsion mesuré en minimisant un \ 2 . La figure ( |3Ï| ) donne un exemple d'une 
telle reconstruction. 

8.3.5 Variance et skewness de la convergence locale 

Dans un contexte cosmologique la convergence locale est donnée par l'ensemble 
des concentrations de masse sur les lignes de visée. Ainsi k dans la direction 7 
est donnée par (dans le cas d'une lentille mince ou par linéarisation), 

«(7) = ~\^j dz s n(z s ) J àx D{Z ^ iz) S masss (x> 7) (1 + zix)) (315) 

où z s est le redshift des sources, et n(z s ) leur distribution en redshift (normalisée 
à l'unité), x es t la distance radiale et D les distances angulaires qui intervien- 
nent (exprimées en unités c/Hq). En général les distances cosmologiques qui 
interviennent sont différentes, sauf dans le cas d'une section spatiale plane. 

On voit que la convergence locale est formellement très similaire à la den- 
sité projetée dans un catalogue angulaire. La différence est que la fonction de 
sélection n'est pas ici normalisée. En utilisant l'approximation des petits angles 
on peut calculer les premiers moments de cette distribution. 
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Figure 32: La variance (panneaux de gauche) et la skewness (panneaux de 
droite) de la convergence locale. En haut est donnée la dépendance avec l'échelle 
angulaire pour deux populations différentes de sources (z s = 1, lignes continues), 
(z s = 2, lignes pointillées). Au centre est représentée la dépendance avec ilo et 
en bas la dépendance conjointe avec Qq et Ao pour une population de sources à 
z s = 1 (figures extraites de Bernardeau el al. 1997). 
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On trouve que (Bernardeau, van Waerbeke & Mellier 1997) 
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(316) 



et 




2 = -40 fi, 



'0 



0.8 



-1.35 



(317) 



Il est sans doute utile de faire un certain nombre de commentaires, 

• La variance dépend non seulement de l'amplitude des fluctuations (car- 
actérisée par as) mais aussi de fio- 

• La skewness est indépendante de l'amplitude des fluctuations (comme on 
s'y attend pour des conditions initiales Gaussiennes), mais dépend de Qo. 
La mesure conjointe de la variance et de S3 est nécessaire pour lever la 
dégénérescence entre as et fio . 

• Ces deux quantités dépendent de la distribution en redshift de la popula- 
tion de galaxie sources utilisée pour faire les mesures. 

Pour l'instant aucune mesure de ces quantités n'a été faite. On en est à des 
études de faisabilité et des investigations sur les erreurs systématiques possibles 
(projet DESCART avec Y. Mellier). Parmi les sources d'erreurs systématiques 
possibles on peut noter, 

• les couplages non-linéaires des lentilles entre-elles en cas de lentilles mul- 



• l'existence de couplages observationels supplémentaires dues en particulier 
à l'effet d'amplification qui peut changer la population d'objets sélectionnés; 

• les effets de clustering des galaxies sources. 

Ces sources d'erreur affectent la robustesse de l'interprétation cosmologique des 
cartes produites. A cela il faut encore ajouter une bonne détermination des 
rcdshifts des sources, mais aussi s'assurer que les cartes produites ont bien une 
origine cosmologique, autrement dit que la qualité optique instrumentale, le 
suivi du télescope, etc, permettent bien la réalisation de cartes fiables. C'est à 
ce jour là que les efforts les plus importants sont faits. 

9 Vers le régime fortement non-linéaire 

9.1 Données expérimentales 

Les données expérimentales suggèrent que la hiérarchie obtenue avec la théorie 
des perturbations se prolonge dans le régime non- linéaire. C'est particulièrement 
vrai dans l'espace des redshift (voir la figure B71) . 



tiples; 
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Par exemple on sait que la fonction de corrélation à trois points des galaxies 
peut se décrirent par le produit de paires de fonctions de corrélation à deux 
points (Groth & Peebles 1977), Le., 

£(ri,r 2 ,r 3 ) = Q [£(r X) r 2 )£(r 2 , r 3 ) + sym.] (318) 

où Q est de l'ordre de l'unité. Dans le cas de la théorie des perturbations 
Q est en fait une fonction sans dimension des angles et des distances entre 
points. Cette dépendance semble disparaître au profit de cette forme plus simple 
dans le régime non-linéaire. Des tentatives ont été faites pour mesurer des 
fonctions de corrélation d'ordre plus élevé mais cela devient très ardu à cause 
des dépendances géométriques. 

Existe-t-il une explication naturelle à ces propriétés? 

9.2 Les solutions auto-similaires 

On aborde ici un régime qui est très non-linéaire et une description avec un fluide 
avec un seul flot serait donc inadaptée. Il faut revenir à l'équation d'évolution 
dans l'espace des phases, 

-/(x, P) i)+-^- /(x,p,t)-mV x .$(x)- /(x,p,t) = 0(319) 
at m a z ax dp 

A<S>( X ) = ——J /(x,p,i)d 3 p. 

On peut facilement montrer qu'il existe une famille de solutions possibles dites 
auto-similaires qui ont la dépendance en temps suivante, 

/(x,p,t) = i- 3 ' 3 /(x/t a ,p/^) (320) 

où les deux constantes satisfont, 

/3 = a + l/3. (321) 

Cette relation s'obtient assez facilement en réécrivant l'équation (|319| ) avec le 
changement de variable u = x/t a et v = p/t 13 . C'est la condition pour que 
la dépendance en temps se factorise. Il est impotant de remarquer qu'alors les 
fonctions de corrélations ne sont fonctions que de u: elles ne dépendent plus 
explicitement du temps. 

Je donne ici les relations de similarité pour le cas d'un Univers Einstein-de 
Sitter (Davis & Peebles 1977). Il existe des formes similaires pour le cas d'un 
univers ouvert dans la limite fi <C 1 (Balian & Schaeffer 1988). On a vu que la 
virialisation des objets entraînait un découplage de leur dynamique par rapport 
à l'expansion générale. Cette propriété se retrouve au niveau de la fonction de 
corrélation à deux points. On s'attend en effet à ce que lorsque cette fonction 
devient grande elle tende simplement à s'opposer directement à l'expansion. Le 
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nombre de voisins d'un objet donné doit en effet rester constant dès que l'on se 
trouve au sein d'un objet très compact, en équilibre gravitationnel (comme une 
galaxie..). Dans la limite non-linéaire (&{r) S> 1), la densité au voisinage d'un 
point s'écrit, 

p , (x)=p(t)&(x,t). (322) 

Il faut alors que Ç(x, t) se comporte comme a 3 à une distance réelle constante 
pour compenser les effets de dilution causés par l'expansion, soit, 

&{x,t) ~ a 3 £(ax). (323) 

On peut généraliser cette idée aux fonctions de corrélation de plus grand ordre 
et on obtient naturellement, 

£ p (xi, . . . ,Xp,i) ~ a 3(p_1) i P {ax!, . . . , ax p ). (324) 

Si on combine ces lois de conservation qui traduisent les hypothèses de corrélation 
stable avec la forme générale attendue pour les solutions auto-similaires, on ob- 
tient des lois d'échelle pour ces fonctions de corrélation, 

6(Ax,t) = X-i&fat) (325) 
e p (A Xl ,...,Ax p ,t) = A-^-^^X!,...^,*) (326) 

où 7 est lié à l'exposant a par 

7 = 6/(2 + 3a). (327) 

Il suffit d'écrire que a 3 Ç(ax) — Ç(x/t a ). De plus si on suppose que la solution 
auto-similaire est toujours valable, du régime linéaire au régime non-linéaire, 
(ce qui stricto sensus implique le spectre de puissance initial suit une loi de 
puissance), alors on peut relier l'exposant 7 à n: (Davis & Peebles 1977), 

| = ,-(«+ 3 > t^ => a = ^ => 7 = 9 ~ ± ^- (328) 
^\t a J 3n + 3 ' 5 + n v ; 

Cette relation est finalement assez bien vérifiée et elle traduit une propriété 
importante qui est que la dépendance spatiale de la fonction de corrélation non- 
linéaire doit pouvoir se reconstruire à partir du spectre initial. 

9.3 La transformée linéaire-nonlinéaire à la Hamilton et 
al. 

Cette idée a été exploitée à son maximum par Hamilton et al. (1991) qui ont 
montré qu'il existait une relation, au moins à peu près universelle, qui reliait la 
fonction de corrélation linéaire à la fonction de corrélation non-linéaire. Pour 
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écrire cette relation il faut écrire une espèce de relation de conservation de la 
matière, 

r 3 (1 + £(r, t)) = constante (329) 
alors £ (r, t) est une fonction de Çnneaire au point r/ de l'espace réel auquel on a 

rf = r 3 (l + £(r,t)). (330) 

Dans le régime linéaire, quand £ est petit, r; et r sont évidemment confondus, 
mais ce n'est plus le cas par la suite. L'idée vérifiée par Hamilton et al. est qu'il 
existe une fonction H (une simple fonction pas une fonctionnelle) telle que, 



£(r,t)=W 



slincaircV h ûqnitialj 



a 2 (iinitial) 

Numériquement ces auteurs ont obtenus, 



(331) 



= z + 0.358 z 3 + 0.0236 x 5 
{ ' 1 + 0.0134 x 3 + 0.00202 x 9 / 2 ' 1 ' 

Cette idée a été plus récemment abondamment reprise pour obtenir des trans- 
formées plus précises, et valable pour toutes les cosmologies, s'appliquant à P(k) 
(Peacock & Dodds 1994, Jain et al. 1995). 



9.4 Les modèles hiérarchiques 

On a là une description phénoménologique du comportement de la fonction à 
deux points (ou de manière équivalente du spectre de puissance), mais évidemment 
cela ne constitue pas une description complète du champ non-linéaire. Je re- 
viens alors aux propriétés d'échelle trouvées pour les fonctions de corrélation 
des grands ordres. 



9.4.1 Décomposition en arbres des fonctions de corrélation 

Un moyen de les réaliser 'à la main' est de supposer que la fonction de corrélation 
à l'ordre p s'écrit comme des produits de p — 1 fonctions de corrélation à deux 
points. C'est d'ailleurs, comme on l'a vu, un moyen assez juste de décrire la 
fonction de corrélation à 3 points. 
Cela revient à dire que, 

), (333) 

config. a liens, (ij) 

où les configurations sur lesquelles cette somme est faite correspondent à toutes 
les manières possibles de connecter p points ensemble avec des diagrammes en 
arbre, les liens sont l'ensemble des liens qui apparaissent entre 2 points (i,j). 
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9.4.2 Contraintes sur les Q p ^ a 

Un des problèmes auquel on est confronté quand on veut construire explicite- 
ment de tels modèles hiérarchiques est qu'on n'est pas libre de prendre les 
paramètres Q Pt a qu'on veut. Il y a des conditions de positivité qui donnent 
un certains nombre de contraintes. Par exemple il faut que la fluctuation du 
nombre de voisins d'un point donné soit positive. Quand £ est grand le nom- 
bre moyen de voisins dans un volume V est essentiellement, nV ^ si la densité 
de point est n tandis que le nombre moyen carré est donné par la fonction de 
corrélation à trois points moyennée, 

(N v ) « nVl (334) 

{N 2 V ) « n 2 V 2 l 3 . (335) 

Si on écrit que (N%\ — (N v ) est positif on obtient, 

Q 3 > (336) 

On peut généraliser ce genre de contraintes et obtenir des conditions sur l'ensemble 
des Q Pt a (Fry 1984b), mais ce ne sera toujours que des conditions nécessaires et 
non suffisantes. Autrement dit on est jamais sûr que le modèle qu'on choisit est 
réalisable par un processus aléatoire quelconque. Si on part de modèles de frac- 
taies, par exemple celui donné par les points occupés par une marche aléatoire 
avec une loi de Levy, on n'obtient pas un modèle réaliste. 



9.4.3 Les modèles 

A priori les coefficients Q p ^ a peuvent dépendre tout aussi bien de l'ordre p que 
de la géométrie a de l'arbre. Différentes hypothèses ont été explorées dans la 
littérature. Saslaw et Szapudi ont fait de nombreuses analyses en supposant 
que les paramètres Q ne dépendent que de p mais dans la suite je vais plutôt 
suivre les hypothèses de Bernardeau & Schaeffer (1992)! Pour réduire le nom- 
bre de paramètres on peut faire l'hypothèse que les paramètres Q p , a s'écrivent 
comme des produits de vertex i/ p , un facteur v p par vertex qui apparaît dans la 
représentation diagrammatique de la fonction àp-points. Alors toutes des fonc- 
tions de corrélation peuvent être construites à partir de la fonction de corrélation 
à 2 points et de la fonction génératrice des vertex Q. Je la note de cette faon par 
analogie avec les résultats de la théorie des perturbations. La fonction Gs de la 
théorie des perturbations peut guider notre choix (c'est d'ailleurs un choix pas 
si mauvais) mais il n'y a a priori aucune de raison de penser que ces fonctions 
sont identiques. 



9.5 Les probabilités de comptage 

Pour calculer les probabilités de comptage on a besoin de connaître la fonction 
génératrice des cumulants. Il est facile de voir que la relation d'échelle (333) 
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induit des propriétés dites hiérarchiques entre les moyennes géométriques^] des 
fonctions de corrélation, Ç p (R), dans des boîtes de rayon R, 

l p (R) ~ S p f -\R), (337) 

où S p sont des paramètres indépendants de R. On définit alors la fonction (p 
comme la fonction génératrice des paramètres S p de la même manière que pour 
le cas perturbatif, 

^( y )=£-S p LJ^_ . (338) 
P =i 

La relation entre les S n et les v p n'est la même que celle donnée dans le cadre 
des calculs perturbatifs que dans une approximation de type champ moyen (qui 
est très bien vérifiée). L'hypothèse de comportement hiérarchique implique que 
la fonction génératrice tp est indépendante de a du moins quand celui-ci est 
grand (on a vu que la théorie des perturbations permettait de calculer la limite 
a — > de cette quantité). Dans une longue étude Balian & Schaeffer (1989) 
ont examiné les propriétés génériques induites par une telle hypothèse sur les 
propriétés des probabilités de comptage. 

9.5.1 Propriétés génériques de la distribution de densité 

On sait que la fonction de probabilité de densité est reliée à la fonction génératrice 
ip(y) à travers, la transformée de Laplace inverse, 

f ^ ■ *1 d*. (339) 



Dans la suite je ne vais pas regarder les effets discrets, mais supposer qu'on 
travaille dans la limite continue (on peut de toute manière retrouver les expres- 
sions des probabilité de comptage par simple convolution avec une distribution 
de Poisson). La difficulté maintenant est de trouver des contraintes sur ip (y). 
Cette fonction est loin de pouvoir prendre une forme arbitraire: il faut que 
toutes les quantités physiques qu'on peut considérer et calculer avec cette fonc- 
tion aient un sens. 

Un cas particulier important est la probabilité de vide -P(O), probabilité de 
ne trouver aucun objet dans un volume de rayon R si le champ de densité est 
décrit par des particules de densité n. Cette probabilité est donnée par 

P(0) = Jdôp(ô) exp[-nV(l + ô)] , (340) 

V = ^R 3 . (341) 



6 les moyennes géométriques des fonctions de corrélation s'identifient naturellement avec les 
cumulants de la distribution de densité filtrée dans cette même boîte 
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En utilisant ( |339| ) on obtient une expression qui s'intègre directement sur S et 
qui donne, 



P(0) 



1 



exp 



<p(y) 



(342) 



„ 27ri y — n V a 2 

Finalement, P(0) est donné par la valeur du résidu et s'exprime simplement en 
fonction de la fonction génératrice ip(y) (White 1979, Balian & Schaefïer 1989) 



P(0) = exp -nV 



y 



= nVa 2 . 



(343) 



Cette expression^] montre en tout cas que (p(y)/y doit être fini et positif pour 
des valeurs arbitrairement grandes de y. Il est alors naturel de supposer que, 



<p(y) 



ay 1 -" »>1, 



(344) 



avec < lu < 1. Les valeurs de y positif correspondent à une série alternée pour 
tp. Par contre quand y est négatif tous les termes de la série ont le même signe 
et il n'y a pas de raison que <p(y) converge toujours. On note y s la plus petite 
(en valeur absolue) valeur de y où <p(y) est singulier, et on peut paramétriser la 
forme de (p(y) par, 

p(y) - p. ~ a. (y - (345) 

Avec ces hypothèses on peut entreprendre le calcul de p(ô). Il apparaît deux 
régimes, l'un correspondant aux régions plutôt sous-denses, l'autre aux régions 
sur-denses. Ici je vais m'intéresser exclusivement au deuxième. Il s'obtient on 
supposant que y ne prend que des valeurs finies ou petites ce qui permet de 
développer une partie des termes sous l'exponentiel et d'avoir, 

dy <p(y) 



P (ô)dô 



27ricr 4 cr 2 



■ exp 



yô 



dô. 



(346) 



On voit que la fonction de distribution peut s'exprimer en fonction d'une unique 
intégrale, h(x), 



p(ô)dô 
h(x) 



dô 



a' 

+ioo 



ày 
2ni 



tp(y)exp(yx). 



Les propriétés de <p(y) induisent les propriétés suivantes pour h[x), 



h{x) 
h(x) 



1 



r( w ) 

a s 
T(-u)s) 



x < 1 

_1 exp(y s x), x > 1 



(347) 
(348) 

(349) 
(350) 



7 On retrouve le cas Poissonnien quand ip(y) = y, c'est à dire quand on a annulé toutes les 
fonctions de corrélation 
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Log 10 O) 

Figure 33: Evolution de la forme de la fonction de distribution de probabilité 
de la densité locale en fonction de la variance. Les cas a = 1, 10 et 100 ont été 
tracés pour des cumulants correspondant à n = — 2. 



On voit donc que la fonction h{x) (et par là même la fonction p(S)) est constru- 
ite par assemblage d'une loi de puissance quand x est petit, et d'une coupure 
exponentielle à grand x (voir figure |33| ). 

Il y a bien sûr une part d'arbitraire dans ce résultat (existence et forme de 
la singularité pour ip(y), ...), aussi il est intéressant de voir ce qui peut être 
construit de manière générique. 

On a vu que dans le cas de la théorie des perturbations, ip(y) pouvait être 
construite à partir d'une fonction génératrice de vertex, G(j). Le comportement 
assymptotique de <p(y) est lié au comportement asymptotique de Q{t) pour r 
positif. Plus spécifiquement si 

g(T)~cT- K (351) 

alors 



" = ^> (352) 
(l-uj)a ^k-^c 1 -^. (353) 

Par exemple pour Q obtenu à partir de l'effondrement sphérique on a k = 3/2, 
ce qui donne u — 3/7 « 0.43. 

L'hypothèse d'une composition en arbres se révèle très féconde pour la partie 
singulière de f(y). En effet l'équation implicite en t induit très naturellement 
une singularité en (p sans faire d'hypothèse particulière sur £7(t). Ce point 
singulier est donné par la solution du système, 

V> =-3*ky (354) 

T s = ç"\ts) • (355) 
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Figure 34: A gauche, détermination de l'indice effectif n e ff à partir des premiers 
moments S p , p — 3,4,5 pour nn n = —1. La courbe pointillée correspond au 
fit proposé dans le texte. A droite les probabilités de comptage mesurées P(N) 
sont comparées aux prédictions théoriques obtenues avec ce n e s- Les simulations 
ont été examinées à différents temps et avec différents rayons de filtrage. 



Alors on a 

<P<V) -<P.~ §(» Vsf /2 Q'(r s )Ç"(r s ) ( 2g, ^) (Ts) ) V2 • (356) 

Dans le cas du modèle de théorie des perturbations cette singularité semble en 
effet jouer un rôle actif. On a vu dans la section précédente quelles étaient les 
valeurs décrivant cette singularité en fonction de l'indice spectral. 

9.5.2 La théorie des perturbations 'étendue' 

Quel modèle pour f(y)7 Au delà du régime quasi- linéaire une idée est de suivre 
au plus près le modèle suggéré par la théorie des perturbations. Formellement 
on peut voir la forme fonctionnelle de tp(y) comme étant paramétrée par l'indice 
spectral n. On peut donc 'oublier' la signification physique de ce paramètre et 
le supposer libre pour voir s'il est possible de l'ajuster pour décrire les résultats 
numériques, aussi bien au niveau des paramètres S p que de la forme des P(N). 
Ce travail a été mené à bien par Colombi et al (1997). L'idée de départ est 
la suivante: les mesures des premiers S p fournissent une estimation de n c g. en 
supposant que 

S p (n cB .) = S~ c . (357) 

Le résultat remarquable de Colombi et al. est que les n e ff. ainsi mesurés ne 
dépendent pas (ou peu) de l'ordre p (pour un spectre initial de forme donné). 
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Table 2: Paramètres utilisés dans le fit ( 358 ). 



^lincar 


"nonlinear 


nonlinear 


nonlinear 


oo 


T 


-2 


-9.5 


-12.4 


-7.22 


1.6 


1.4 


-1 


-3 


-3.8 


-2.24 


1.4 


1.2 





-1.2 


-1.6 


-0.86 


1.25 


0.6 


+ 1 


-0.85 


-1.17 


-0.57 


0.7 


0.3 



On le voit sur la partie gauche de la figure (p4|). Du coup les probabilités de 
comptages qu'on obtient en utilisant le n e g ainsi mesuré reproduisent très bien 
celles qui sont obtenues numériquement (partie droite de la figure [Ï4j) . 
Le paramètre n e g est donné par 

n e s = nun + (n non _ lin - n lin ) — — — — (358) 

x T + x T 

x = exp[log 10 (cr 2 /cro)]. 

Cette transition est décrite par la valeur de n non -iin atteinte pour les grandes 
valeurs de er, le moment de la transition aç, et la largeur de cette transition, r. 
Les paramètres trouvés sont donnés dans la table (0). 

Les figures (|34|) montrent que la descriptions des moments et des P{N) 
en terme d'un simple indice effectif sont effectivement prometteuses. Il n'y a 
cependant aucune justification théorique en faveur d'un tel comportement. 



9.6 Distribution de matière et distribution de lumière: 
vers une théorie unifiée 

Ultimement ce qu'on veut c'est une description complète des objets présents 
dans l'Univers, en particulier les propriétés de corrélation des puits de po- 
tentiel qui se sont créés. Ils sont implicitement contenus dans la description 
des fonctions de corrélation, aussi compliquées qu'aient pu être leurs histoires 
(accrétion, merging, ...). Potentiellement on doit donc faire beaucoup mieux que 
des théories à la Press et Schechter (1977) qui font des hypothèses drastiques 
sur la dynamique. 

La mise en œuvre de ce calcul est cependant assez complexe. L'identification 
mathématique d'un puits de potentiel est une chose qui est loin d'être aisée. 
Cependant on est maintenant dans une situation assez favorable: les objets 
que l'on recherche ont des contrastes en densité très grands et leur détection 
ne doit dépendre que très peu du critère de sélection utilisé. Et en fait on va 
supposer que la simple probabilité d'avoir localement un contraste en densité 
suffisamment grand va suffire à définir un objet. 

En conséquence de quoi la fonction h(x) donne une description de la forme de 
la distribution de masse d'objets de taille donnée. Les propriétés de corrélation 
de ces objets seront données par les probabilités de comptage conjointes (à deux 
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points si on veut la fonction de corrélation à deux points, à 3 points si on veut 
celle à 3 points, etc.). 

L'hypothèse d'une décomposition en arbres est alors encore extrêmement 
féconde pour faire le calcul de ces probabilités conjointes. 

9.6.1 Le biais 'dynamique' 

Il faut commencer par calculer la distribution conjointe, p(ôi,Ô2), probabilité 
d'avoir la densité ô\ et la densité S 2 conjointement en deux points aux posi- 
tions xi et x 2 . Pour obtenir une telle distribution il faut calculer la fonction 
génératrice ip(yi,y 2 ) des cumulants conjoints pris en deux points distincts. Plus 
précisément on peut définir C pq par 

<^( Xl ) 5« (x 2 )) c = C pq (5l) p+q - 2 (Mxi) Mx 2 )). (359) 

Cette relation définit des paramètres finis C pq dans la limite où la distance 
|xi — X2I est beaucoup plus grande que l'échelle de filtrage R. A partir de là on 
définit la fonction génératrice double, 

<p(vi,V2) = 2^2^ c m — j ~\ — • ( 36 °) 

p=0 q=0 F ' H ' 

en ayant posé, 

C p o = S p , Coo = 0. (361) 

La décomposition exacte en arbres implique que les facteurs C pq ont en plus 
une propriété de factorasibilité, c'est à dire que, 

C pq = C p C q pour p>l et q > 1. (362) 

C'est une propriété assez facile à voir d'un point de vue diagrammatique: le 
poids qu'on mets d'un côté du diagramme ne dépend pas de la topologie qu'on 
met de l'autre coté. De plus, compte tenue de la remarque qui a été faite plus 
tôt, la fonction r(y) s'identifie à la fonction génératrice des vertex à une patte 
externe. Plus précisément on a, 

r(2/ ) = _]rc p ^— , (363) 
p=i 

quand r(y) est solution du système donnant <p(y). Cette relation fonctionclle 
est extrêmement intéressante mais dépend de manière explicite de l'hypothèse 
de décomposition en arbres des fonctions de corrélation. 

Si on exploite cette propriété on peut calculer la probabilité conjointe p{5\ , ô 2 ) 
qui s'écrit, 

p(ôi,ô 2 ) = ^h(x 1 )h(x 2 ) [l + 6(x 1 )6(x 2 )6(xi,x 2 )] (364) 
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avec 

/ioo 
— r(y) exp(xy). (365) 
-ioo 27Tl 

Autrement dit les régions les plus denses continuent à suivre une loi hiérarchique. 
Le facteur b(x) s'identifie au facteur de biais entre les pics. Cela implique que le 
biais est a priori une fonction de x = uniquement (et pas de 5 par exemple). 
La forme de b(x) n'a pas a priori une forme indépendante des paramètres des 
modèles. Cependant la singularité qui apparaît naturellement pour tp(y) induit 
un comportement spécifique pour b(x), 

b(x) ~ x quand x » 1. (366) 

Le résultat générique est donc que le biais croit avec la rareté. Ce résultat 
obtenu à partir d'une description aussi réaliste que possible du régime non- 
linéaire prolonge d'une certaine manière les calculs de Bardeen et al. (1986) 
menés dans un champ Gaussien. La principale conséquence de ces calculs est 
qu'a priori les puits de potentiel qui se forment sont effectivement biaisés. 

9.6.2 La hiérarchie des corrélations pour les pics de densité 

Le paragraphe précédent présente les résultats pour la fonction de probabilité 
à deux points. On peut en fait généraliser ces résultats pour des fonctions 
de corrélation d'ordre donné. Dans l'hypothèse de fonctions de corrélation en 
arbre, on peut définir la fonction génératrice des arbres à un nombre arbitraire 
de pattes. Cela permet de calculer l'ensemble des propriétés de corrélation 
jusqu'à un ordre donné. Les résultats sont les suivants, 

• Les propriétés hiérarchiques des fonctions de corrélation aux grands ordres 
sont conservées. 

• Les paramètres S p correspondant sont des fonctions de x uniquement. 

• Dans la limite x — > oo, ces paramètres ont une limite finie donnée par 

S p {x)^>p p - 2 quand x -> oo. (367) 

Ces résultats sont assez généraux et dépendent finalement peu des hypothèses 
précises sur les fonctions de corrélation. Ils supposent des propriétés de régularité 
de la fonction génératrice des vertex. Une conséquence, peut-être importante, 
de ces résultats est que les S p des pics rares devraient être les mêmes qu'on soit 
dans le régime quasi-linéaire ou dans le régime fortement non-linéaire. C'est 
peut-être le début d'une explication pour les observations qui sont faites. 

Ces résultats montrent que même en l'absence de mécanismes non-gravitationnels 
il est possible d'avoir des effets de biais. Il faudrait, pour faire mieux, connaître 
les solutions explicites du régime non-linéaire. 
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10 Conclusions 



Dans ce petit tour d'horizon de la cosmologie j'ai été loin d'être exhaustif. Dans 
le domaine de l'évolution des grandes structures je n'ai pas mentionné bon 
nombre de domaines de recherche en plein développement, comme 

• L'étude des flots cosmiques à grande échelle, et leur utilisation pour mesurer 
les paramètres cosmologiques; 

• L'étude des amas de galaxies et de leur contenu aussi bien en matière 
noire (études dynamique, reconstruction de masse par effets de lentilles 
gravitationnelles..), qu'en matière baryonique (galaxies, rayonnement X..); 

• L'étude des absorbants et des objets de grand z; 



Du point de vue du physicien théoricien un domaine de recherche de prédilection 
est inconstablement la cosmologie primordiale, lieu de rencontre entre entre la 
physique des hautes énergies et la cosmologie observationnelle. 

Cependant la formation des grandes structures recelle un certain nombre de 
problèmes ouverts digne d'intérêt, 

• La théorie des perturbations est loin d'avoir livré tous ses secrets. En 
particulier il serait extrêmement utile de comprendre les corrections en 
boucles. 

• De nombreux aspects du régime non-linéaire ne sont pas compris. Par 
exemple on ne sait pas décrire la transition vers le régime multiflots qui 
conduit à la virialisation de la matière dans les puits de potentiel. 11 
serait aussi très intéressant de pouvoir exhiber une solution explicite des 
équations dynamiques dans le régime nonlinéaire (même si ce n'est qu'une 
forme assymptotique) . 

• L'exploitation des catalogues tridimensionnel ou bidimensionnel notam- 
ment avec la mise en évidence de propriétés non-Gaussiennes n'est pas 
encore optimale. La question se pose par exemple pour les cartes de dis- 
torsion gravitationnelle. 

• La relation entre baryons et matière noire est un domaine qui n'a pra- 
tiquement pas été exploré. On ne dispose pour l'instant que d'expériences 
numériques pour essayer de comprendre ce qui se passe. 
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